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Abstract 



The chemotaxis pathway in the bacterium Escherichia coli allows cells to detect changes in external 
ligand concentration (e.g. nutrients). The pathway regulates the flagellated rotary motors and hence the 
cells' swimming behaviour, steering them towards more favourable environments. While the molecular 
components are well characterised, the motor behaviour measured by tethered cell experiments has 
been difficult to interpret. Here, we study the effects of sensing and signalling noise on the motor 
behaviour. Specifically, we consider fluctuations stemming from ligand concentration, receptor switching 
between their signalling states, adaptation, modification of proteins by phosphorylation, and motor 
switching between its two rotational states. We develop a model which includes all signalling steps in 
the pathway, and discuss a simplified version, which captures the essential features of the full model. 
We find that the noise characteristics of the motor contain signatures from all these processes, albeit 
with varying magnitudes. This allows us to address how cell-to-cell variation affects motor behaviour 
and the question of optimal pathway design. A similar comprehensive analysis can be applied to other 
two-component signalling pathways. 
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Introduction 

Biological systems sense stimuli from their environ- 
ment using cell-surface receptors, and process this 
information to make reliable decisions, e.g. where 
to move, or whether to divide or to express new en- 
zymes. Typically, intracellular signalling molecules 
are activated by modification, e.g. phosphorylation 
and methylation, and interact in complicated bio- 
chemical reaction networks. Biochemical reactions 
rely on probabilistic collisions of a limited num- 
ber of molecules. Hence, the number of signalling 
molecules fluctuates with time, making signal pro- 
cessing noisy. The abundance of noise sources in a 
cell is in stark contrast to the remarkable accuracy 
with which cells are known to respond to minute 
amounts of chemical concentr ation, including grow- 
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ing axons and immun e cells (jSvkulev et al 
iMortimer etafl , |2010|) . 

The high biological relevance of noise has 
widely been recognised and studied exten- 
sively in gene expression (Elowitz et al. . 2002 



Pedraza and van OudenaardenL 20051 IPaulsson . 
20051: lAcar et all l2010t lEldar and Elowitzl |201qF . 



In contrast, noise in signal transduction is not well 
characterised, despite its importance for accurate 
sensing and cell-decision making. Examples of eu- 
karyotic systems, in which signalling noise has been 
considered include the ultrasensitive thre s holdin g 
cascades ( Thattai and van Oudenaarden j_ 20021) 



pheromone sensing i n yea st (jColman-Lerner et al 



20051 iTavlor et all 120091) . signal tr ansduction in 



photoreceptors ( Detwiler et ah L 2000h and feedback 
loops for noise suppress ion ()Hornung and Barkail 



2008 ; Lestas et al. . 2010h . Furthermore, signalling- 
noise has been considered in part s of bacterial 



pathways dEmonet and Cluzell. 2008 : Mehta et all 
20081: iMora and Wingreenl . l2010h~ However, analy- 



ses have either been not comprehensive, or signal 
and noise transmission have not been compared ex- 
plicitly. 

An important class of signalling pathways are 
the bacterial two-component systems, including 
hundreds of pathways responsible for wide rang- 
ing functions such as sensing of and responding to 
nutrients, os molarity, antibiotic s , as w ell as quo- 
rum signals (jLaub and Goulianl 120071) . A partic- 



ularly well characterised example is the chemo- 
taxis pathway in E. coli (Fig. [1]), allowing cells 
to swim towards nutrients and away from toxins 
with high sensi tivity over a wide range of ambient 
concentrations (Berg, 2000; Falkc and Hazclbaucr, 



Baker et al.l 



20011: ISouriikL 120041 iWadhams and Armitagd . 12004 ; 

20061 ). Specifically, the kinase 



CheA autophosphorylates when receptors are ac- 
tive and passes on phosphoryl groups to the re- 
sponse regulators CheY and CheB. Phosphorylated 
ChcY (CheY-P) modulates the probability of coun- 
terclockwise (CCW) or clockwise (CW) rotation of 
the motor. The rotational directions of motors cor- 
respond to the two swimming modes of the bac- 
terium, namely smooth swimming and tumbling, 
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Figure 1: Chemotaxis pathways of E. coli. The 
pathway consists of transmembrane chemorecep- 
tors, which constantly undergo molecular transi- 
tions between their on (black) and off (white) con- 
formations. Receptors signal to CheA, which sub- 
sequently autophosphorylates. The average re- 
ceptor activity is modulated by ligand binding 
and unbinding, as well receptors methylation and 
demethylation. CheA phosphorylates the response 
regulator CheY, which defuses through the cell 
and binds to the rotary motors. Upon CheY-P 
binding, motors switch from their default state of 
CCW rotation (i.e. running mode) to CW rota- 
tion (i.e. tumbling mode). In addition, CheY-P 
is dephosphorylated by its phosphatase CheZ. Re- 
ceptor methylation is catalysed by CheR, which 
preferentially modifies inactive receptors. Recep- 
tor demethylation is catalysed by CheB, which is 
activated by phosphorylation, and modifies prefer- 
entially active receptors. 



respectively. Adaptation, i.e. the reversal of the 
effect of changes in the ligand concentration, is 
mediated by reversible receptor methylation and 
demethylation, catalysed by enzymes CheR and 
phosphorylated CheB (CheB-P), respectively. 

Using the E. coli chemotaxis pathway as an ex- 
ample, we are interested in the behaviour of the ro- 
tary motor, i.e. the cell's final output, and how its 
rotation is affected by signalling and noise. Specif- 
ically, there are several fundamental questions we 
would like to address: 

Firstly, what types of signals are transmit- 
ted and what types are attenuated by the path- 
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way? Early work showed that the system re- 
spon ds to the t i me-de rivative of the input sig- 
nal (|Block et al.l . 119821) . A number of research 



groups have measured the av eraged response of 
cells to chemotactic signals dBlock et all Il982 



Seeall et all Il986t IShimizu et alll2010l) . and found 



that slowly, as well as rapidly changing input sig- 
nals are not transmitted by the pathway. The 
response to slowly changing signals is attenuated 
by adaptation, which reverses the activation by 



ligand binding (Block et all 119831 iTu et all 2008 



[Sh imizu et ail l2010h . Rapidly changing signals 



were conjectured to be attenuated by a thi rd-order 
filter dBlock et allll982HSegall et allll98fih . While 
the phosphorylation dynamics of CheY-P has been 



shown to contribute a first-order filter (Tuetal 



2008), the exact filtering dynamics of the full path- 



way has not been addressed. 

Secondly, how is noise generated, amplified 
or filtered in the signalling pathway, and how 
do different sources of noise affect the motor 
behaviour? The power spectrum, which captures 
the correlations between fluctuations in motor 
behaviour at different time points, was measured 
for wild-type cells and mutant cells lacking the 



chemotaxis signalling pathway ([Korobkova et al 
2004). The spectrum was found to have a large 
low-frequency component in the wild-type cells, 
indicating that there is a dominant noise source 
in the signalling pathw ay wit h long correlations. 
Corobkova et al. ( 2004h and Emonet and Cluzel 



2008]) showed, using simulations of the signalling 
pathway, that the adaptation dynamics plays 
an important role in generating long correla- 
tions. However, they only analysed the signalling 
pathway up to CheY-P. Other studies include 
stochastic simulations of the noisy biochemical 
react ions of the pathway ([Morton-Firth and Brav . 
1998), and addressed the mechanism of mo- 



2006: Meacci and Tu 



tor r ot a tion (|Xing et al 
20091: Mora. Yu. Sowa and Wingreen , 



2009; 



Mor a. Yu and Wingreenl . I2009t Ivan Albada et al 
2009h . inc luding the thermodyna mics of motor 



switching dScharf et all Il998t iTurner et all Il999t 



Tu and Grinstein . 20051) . However, noise gener- 



ation, filtering and amplification has not been 
addressed systematically for the various levels of 
the signalling pathway from chemoreceptors to 
motors. 

Finally, how reliably are signals transmitted in 
the presence of noise? An important task for the 
cell is to generate an appropriate response to input 
signals in the presence of fluctuations in the input, 
as well as due to noise in the biochemical signalling 
pathway. Furthermore, cell-to-cell variation in pro- 
tein expression influences signal transmission and 
noise filtering. Comparing these two aspects of the 
pathway dynamics, the signal-to-noise ratio, is a 
novel perspective in our present study. 

In the following, we present a mathematical 
model for the chemotaxis signalling pathway. A 
simplified pathway is given in the main text, while 



details of the full pathway are provided in the 
Supplementary Information (ST). We discuss the 
average (deterministic) response of the signalling 
pathway to concentration signals. We analyse the 
noise sources in the signalling pathway and their 
effects. Finally, we vary pathway parameters and 
study how they affect signal and noise transmis- 
sion. We also discuss briefly how our approach can 
be applied to other two-component systems and sig- 
nalling pathways. Introductions to our modelling 
approach (Box 1), the mathematical characterisa- 
tion of signal and noise propagation (Box 2), as 
well as a comparison of E. coli's chemotaxis path- 
way and other two-component pathways (Box 3) 
are also given. 



Results 

Experimental measurements of re- 
sponse and noise spectrum 

The signal propagation in the chemotaxis pathway 
has been characterised by the response to small 
concentration signals (linear response function; see 
Box 2). Specifically, the response has been mea- 
sured at the level of CheY-P using fl uorescence res- 
onance energy transfer (FRET) by [Sh imizu et al 
(|2010h . In that study the system was stimulated 



by a periodic variation of the concentration of at- 
tractant MeAsp. Using a series of frequencies of the 
stimulation, the magnitude (modulus) and phase, 
i.e. the lag between signal and response, of the re- 
sponse was determined. In cell-tether experiments 
of motor rotation, the response to short impulses 
of attr actants was measur ed at the level of t he mo - 
tor bv iBlock et "all (|l982l) and ISeeall et all (| 19861) . 



Such data determines the linear response function 
up to a constant factor. Experimental results are 
shown in Fig. [5] Noise propagation in wild-type 
cel ls has only been s t udied at the level of the motor 
bv lKorobkova et al.l (|2004h . We use the experimen- 
tal response functions to calibrate our model, and 
subsequently study the noise power spectrum and 
signal-to-noise ratio. 



Simplified model for the pathway 

Here we consider a simplified pathway to gain intu- 
ition of the key processes involved. The simplified 
pathway consists of chemoreceptor signalling in re- 
sponse to ligand binding and receptor methylation, 
as well as the rotary motor. Specifically, we use 
stochas tic differential equat ions in a Langevin ap- 
proach (jvan Kampenl . 120071 : see Box 1) to describe 
the dynamics of each type of signalling protein. We 
assume throughout that fluctuations in concentra- 
tion are small, allowing us to describe the average 
behaviour of a signalling molecule by a determinis- 
tic dynamics and the fluctuations around the mean 
by additive noise. 

We assume N receptors form cooperative sig- 
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Box 1 Langevin approach for signalling in the pathway. 



- L R. 



In this paper, we use the Langevin approach to de- 
scribe the noisy dynamics of signalling. This approach 
is based on the assumption that on average the dy- 
namics of signalling follows a deterministic equation, 
an ordinary differential equation (ODE) for the rate 
of change. Consider for instance the following (lin- 
ear) ODE for the rate of change of R(t), which could 
describe the phosphorylated fraction of a molecular 
species, 

dR _ 
dt ~ 

This equation describes how 7? changes due to con- 
stant, basal production with rate k and decays with rate 
t~ 1 R, i.e. proportionally to the number of molecules 
R. The latter term results in the exponential decay of 
R (similar to radioactive decay) with a characteristic 
time constant r (or characteristic frequency r~ 1 ). 

The Langevin approach assumes that fluctuations 
around the average dynamics are small and well de- 
scribed by a noise term r)(i) added to the above (deter- 
ministic) dynamics. The resulting equation is a stochas- 
tic differential equation (SDE), 



— - k 
dt ~ 



l R+T)(t). 



The noise term r/(t) is typically assumed to be Gaus- 
sian and white. The first property means that fluctua- 
tions around the average value of R are assumed to be 
normally distributed. The second property means that 
the noise contains all frequency components or colours, 
and specifically that the autocorrelation function (see 
Box 2) is {r}{t)r]{t')) = Q8(t - t'). Here, Q is the in- 
tensity of the noise and 5(x) the Dirac delta function. 
The angular brackets indicate averaging over numer- 
ous time series. This autocorrelation function indicates 
that fluctuations in the rate at any two different time 
points are not correlated. 

Linearisation of the equations. In order to 
characterise the response to input signals, as well as 
fluctuations in Box 2, we need to linearise and Fourier 
transform the dynamical equations. Here, we illustrate 
this procedure for our simple example. Considering our 
ODE from above, we note that the variable R{t) as- 
sumes a constant steady-state value R* when its rate 
of change is zero, i.e. dR/dt = 0. We obtain the dynam- 
ics of deviations of R(t) from its steady state, e.g. due 
to an input signal, by linearising R(t) = R* + AR(t) 
around the steady state. Inserting this linearisation 
into the above example yields 



d[R* + AR(t)] 
dt 



= k-T~ 1 {R* + AR(t)]. 



As the rate of change of the steady state is zero by defi- 
nition, we obtain for the dynamics of deviations AR(t) 



d[AR(t)\ 
dt 



= -T _1 Ai?(i). 



Considering the SDE, we are interested in the dy- 
namics of fluctuations in R(t), i.e. deviations from its 
average value (R). Hence, we write R(t) — (R) +SR(t), 
with SR(t) fluctuations in R. Inserting into the SDE 
yields 



d[{R) + SR(t)] 
dt 



= k-T- 1 [(R)+5R(t)]+r 1 (t). 



The dynamics of the average (R) is described by the 
deterministic dynamics described by the ODE, whereas 
the dynamics of fluctuations is 



d[SR(t)} 
dt 



^SR^ + riit). 



Fourier transformed equations. Analysis of the 
dynamics is simplified by considering the Fourier trans- 
formed equations, defined for AR(t) as 



AR(cj) 



AR(t)e mt dt. 



Inserting the inverse Fourier transforms AR(t) = 
1/(2tt) AR(uj)e- iujt duj into the ODE for AR{t) re- 
sults in the equality 



-iu>AR(uj) 



_1 Aii(cj) 



i.e. the time-dependent differential equation has been 
transformed into an algebraic equations, which can eas- 
ily be solved for AR(uj). A similar equation can be 
obtained for the SDE for the fluctuations. 



Determining the noise intensity. To discuss 
this aspect, consider the following general SDE: 



dR 
~~dt 



n - r 2 + r)(t). 



The rates r\ and r 2 typically depend on the concen- 
trations of proteins in the signalling network. The 
noise term r/(t) is composed of two terms r)\(t) and 
T72 (t), which are associated with the rates ri and r 2 , 
respectively. We assume rji and r\ 2 to be indepen- 
dent, i.e. {rji(t)ri2(t')) — 0. In general, this is jus- 
tified as different reactions are catalysed by different 
proteins. Using (rij(t)r]j(t')) = QjS(t — t'), the noise 
intensities can be calculated if we make the assump- 
tion that fluctuations are due to so-called birth and 
death processes, i.e. creation and destruction of the 
molecules with average rates r\ and r 2 . Then the 
ass ociated noise intensities are Q i = r \ and Q2 = 
r 2 (|Thattai and van OudenaardenL l2002h . The inten- 
sity of the total noise r](t) is the sum Q = Qi + Q 2 due 
to the independence of the two noises. As forward and 
backward rate are equal at steady state, Q is twice the 
reaction rate in one direction at steady-state. 
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Concentration [Y**] [nM] 



Figure 2: Calibration of the model. (A and B) Fourier transform ed linear response function. (A) 
Magnitude of the response function of the rotary motor measured by iBlock et all {l982) (circles) and 
Segall et al. (|l986f ) (plus symbols). The fit of our model is shown by the solid line. The dashed line 
represents a 3rd-order filter for comparison. (B) Magnitude (left ) and phase (r i ght) o f the response 



20ld ). The symbols 



function at the level of the response regulator CheY measured by (Shi mizu et al. 
are measurements at 22° (circles) and 32° (plus symbols). Grey lines represent the fit of our model to the 
magnitude of the response. (C) Switching rates of the motor from CCW to CW rotation k + (squares) 
and from CW to CCW rotation fc_ (circles) as a fu nction of the concentration of signalling-active mutant 



Y**. A fit using the model of iTurner et al. ([1999) is shown as well (solid and dashed lines; cf. SI). 



nailing complexes, which can switch between an 
active (on) and an inactive (off) state. Their 
activity A is described by th e Monod-Wyman- 
Chan g eux (MWC) mod e l dSouriik and Bere 



given by 



2004 iMello and Tul. [20o5 iKevmer et al 



Endres and Wingreenl. 20061 Endres et al 



2006 



Clausznitzer et alT 2010ll . The activity depend; 



2008 



s on 

the external ligand concentration c at the receptor 
complex, as well as the methylation level M of the 
complex as detailed in Materials and Methods. 

We consider Nc receptor complexes in a cell, and 
assume that each complex signals independently of 
the others. The total activity A c of all receptors in 
a cell is determined by the sum over all signalling 
complexes j. The dynamics of the total activity is 



^ dA dM j 



dA r 



dt ^ dM dt 

l — l 



dA dcj , , 



i.e. the dynamics of the complex activity is af- 
fected by changes in the receptor complex methyla- 
tion level (first term), changes in ligand concentra- 
tion (second term), as well as fluctuations due to 
the switching of the complex between its states (last 
term). All noise terms r\(t) introduced in this sec- 
tion are discussed in Materials and Methods. 

Changes in the concentration originate from 
time- varying input signals (c(t)), as well as fluc- 
tuations due to ligand diffusion. The dynamics of 
the concentration at the jth receptor complex is 



dt 



dim 

dt 



Vc, (t), 



(2) 



where the first term captures average concentration 
changes (indicated by angular brackets (■••)), af- 
fecting all receptors, and the second term describes 
concentration fluctuations at each receptor com- 
plex. 

Adaptation is provided by reversible recep- 
tor methylation and demethylation, whose dy- 
nami cs is described by th e following equa- 
tion ( Clausznitzer et aD . l2010h : 



dMj 
dt 



1r(N - A,) 



lBA] + n M] (t). (3) 



The total methylation level Mj of a receptor com- 
plex j is changed by methylation of receptors in the 
inactive state (first term) and demethylation (sec- 
ond term). This latter rate is assumed to be 
strongly dependent on the receptor complex activ- 
ity as only active receptors are demethylated by 
phosphorylated demethylation enzymes. The last 
term represents fluctuations due to the noisy pro- 
cessivity of the methylation and demethylation en- 
zymes. 

The motor is described as a two-state system 
with CW and CCW rotating states, correspond- 
ing to running and tumbling modes, respectively. 
The dynamics of the probability of tumbling mode 
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(tumble bias) Pew is described by 



dP c 



w 



dt 



k+(A c )(V 



■ Pew )-k-(A c ) Pqw + VPe 



(4) 

where the first term represents the switching from 
CCW to CW with the transition rate k+, the sec- 
ond term represents switching from CW to CCW 
with transition rate fe_ and the third term de- 
scribes temporal fluctuations in switching rates. 
Here, transition rates are modulated by the re- 
ceptor signalling activity A c , whereas in the full 
pathway model CheY-P modulates motor switch- 
ing. These rates have been experimentally mea- 
sured using signalling mutants expressing varying 
amoun ts of constitut i yely a ctive signalling molecule 
CheY (jTurner et all Il999h . The s witching rates 
inclu ding a fit of the model we used ( Turner et all 
19991 : cf. ST) to the data, are shown in Fig. [2] 



Signal propagation 

We consider the response to concentration signals 
at various levels in the signalling pathway to study 
how signals are transmitted to the rotary motor 
(see Box 2 for an introduction to the formal char- 
acterisation of the response). Briefly, an input sig- 
nal Ac(t) is a concentration change relative to a 
constant background concentration cq, affecting all 
receptors equally and representing a "meaningful" 
input to the chemotaxis signalling pathway. Hence, 
the concentration is given by (c(t)) = cq + Ac(t). 
Furthermore, cells are assumed to be adapted to 
the pre-stimulus concentration cq with the various 
levels R of the signalling pathway adapted to their 
steady-state values R* . 

Analytical results for linear response func- 
tions. We can analytically calculate the Fourier 
transformed linear response function from the dy- 
namical equations Eq. (HJ-(H|) without noise (see 
Materials and Methods). Knowing the response 
functions allows us to calculate the response to an 
arbitrary input signal (Box 2). Furthermore, we 
can analyse the filtering of the signal at each level 
of the pathway. The Fourier transformed linear re- 
sponse function for the total activity of all receptors 
in a cell is 



ujm — iuj 



(5) 



The receptor activity is a high-pass filter: The mag- 
nitude of the response function is small for frequen- 
cies lo below ujm = (jr + 3"fBA* 2 )dA/dM, which 
is the characteristic frequency due to adaptation. 
For frequencies above u>m the response function is 
a constant, given by the number of receptor com- 
plexes Nc participating in the response, and their 
sensitivity dA/dc to ligand, evaluated at steady- 
state. The sensitivity is proportional to the recep- 
tor complex size N, i.e. it describes the amplifica- 
tion of the response of a single receptor. 



Similarly, the Fourier transformed response of 
the motor is given by 



xpcwM 



L0 2 



UJ Pc 



(6) 



The motor is a low-pass filter, i.e. its dy- 
namics introduces a frequency-dependent response, 
which is constant below the characteristic fre- 
k*_ of the motor due to 



quency ojp c 



the steady-state switching rates k*L and k*_. The 
parameter u>2 describes the sensitivity of motor 
switching with respect to changes in receptor activ- 
ity (Materials and Methods). At frequencies above 
wp cw the response is reduced. From Eq. © it is 
obvious that receptors and motor are in a cascade: 
The motor response introduces a new filter propor- 
tional to (oji — iuj)^ 1 which simply multiplies the 
response function of the response of the receptor 
activity. The response functions of the full pathway 
including the phosphorylation reactions are shown 
in the SI. 

For further analysis, we can write the Fourier 
transformed linear response function as 



(7) 



where \xr\ is the magnitude and <Pr is the phase of 
the response function, which characterise the am- 
plitude and lag of the response behind the input 
signal, respectively. 

Model calibration. Figure [5] shows experimen- 
tal data for the response function, as well as the 
fits of our full pathway model. For the fit o f our 
model to the data by [Sh imizu et al.l (12010), we 
adjusted only the adaptation rates, as measure- 
ments were restricted to low frequencies. The fit at 
32° C yields the same adaptation parameters as ob- 
taine d from fitting dose-respons e curves of adapting 
cells (|Clausznitzer et all |2010T ) (Fig. H left). The 
adaptation rates for room temperature are one or- 
der of magnitude smaller. Importantly, fitting to 
the magnitude of the Fourier transformed response 
yields a good fit for the phase of the response as 
w ell (Fig.H right). 



Block et alJ (jl982l ) and ISegall et al.l (jl986l) mea- 



sured the response of the motor using impulses 
of attractant. For our fit we adjusted adapta- 
tion and motor switching rate s. Compared to the 
data by [Sh imizu et al.1 (2010) at the same tem- 
perature, adaptation rates are one order of mag- 
nitude larger, i.e. adaptation is faster in these 
experiments. The parameter wp cw of the motor 
switching is 2.1 /s, consistent with switching rates 



of about 1 Hz ( Block et all 11983 ) . It is not clear 



from where the difference in adaptation rates be- 
tween the two sets of experiments originates. Be- 
sides diff erent experimental conditions, it may be 
due to Shimizu et al. (l2010h u sing populat i ons o f 
cells, whereas measurements bv lSegall et all (|l986) 
were done on single cells. Fitted parameters are 
given in the SI. 
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Box 2 Characterisation of signal and noise propagation. 



Signal 

Input | | Output | 



Noise 



Ac 



AR 



Time 



6fl 




Time 



Time 




Frequency 



Frequency 



Frequency 



Signal response. The response of a system to an 
arbitrary small signal is described in terms of a char- 
acteristic function of the system, the linear response 
function \R- The time-dependent response AR(t), i.e. 
the deviation from the adapted state due to a small 
input signal Ac(t), is linear and determined by 

AR(t) = [ X n(t - r)Ac(r)dT. 

Hence, the time course of the response is determined by 
the convolution of the linear response function and the 
input signal. The linear response function describes the 
dynamics of the pathway and the convolution with the 
input signal represents the fact that the current state 
of the syst em is determ ined by the history of the in- 
put signal l|Kubcl.ll957h . The Fourier transform of this 
equation reads more simply 

AR(lj) = xb(w)Ac(w). 

Xr(lo ) is also c a lled the frequency-dependent 
gam (|Detwiler et all . |2000T ). The magnitude 
describes what frequencies of the input signal are trans- 
mitted well, and which ones are attenuated. 

Typically, finite activation rates of the system limit 
the response to rapidly changing input signals, i.e. 
high-frequency signals. In this case, the Fourier trans- 
formed linear response functions falls off at high fre- 
quencies, and the system is called a low-pass filter. If 
low-frequency components of the input signal are fil- 
tered out rather than high-frequencies, the system is 
called a high-pass filter. The chemotaxis pathway is a 
band-pass filter (see figure), filtering out low and high- 
frequency components. 

To obtain a succinct measure for the signalling re- 
sponse due to an input concentration change Ac(i), 
we define AT? 2 the integral over the response over fre- 
quency 



dw|xjj(w)Ac(o;)|' 



Noise. Typically, any system is subject to various 
sources of noise, i.e. random fluctuations in the input, 
as well as from the internal signal processing. This 



is true in particular for biological systems, which rely 
on biochemical reactions and conformational changes of 
signalling molecules, which are intrinsically probabilis- 
tic. 

Fluctuations 8R(t) around the mean value (R(t)) can 
be characterised by their correlations at different time 
points. The autocorrelation function K is defined as 

K(t) = (8R(t)8R(t + r)), 

i.e. the average value the product of fluctuations at two 
time points. It only depends on the interval between 
time points if the dynamics of R(t) is stationary, i.e. the 
mean value {R(t)} and variance {8R 2 (t)} do not vary 
with time t. Averaging over different measurements of 
R(t) is indicated by angular brackets. Typically, corre- 
lations decay with the interval length r between time 
points. Often, the power spectrum is used to charac- 
terise fluctuations. According to the Wiener-Kinchin 
theorem the power spectrum is the Fourier transfo rm 
of the autocorrelation function (|Stratonovichl . IT963I ). 



K(t)e' ut dt. 



For exponentially decaying correlations as in the fig- 
ure, the power spectrum is Lorentzian, i.e. has the 
frequency dependency 



Sr(oj) oc 



1 



LO 2 + U)i 



The power spectrum can be calculated from the abso- 
lute square of the Fourier transform of time series 8R(t) 
measured or simulated over a time interval T, 



Sr(u) = lim 

T—>oo 



(SR(lj)8R*{w)) 
T 



where the Fourier transformation is defined on the fi- 
nite measurement interval T and the average (■} is over 
multiple time series. 

The variance of a stationary process can be calculated 
as the integral of the power spectrum over frequency, 



(8R 2 ) 



_l r 

2W-, 



du Sr(lj) 
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Signal filtering along the pathway. Figure [3] 
shows simulated time courses of the chemotactic re- 
sponse to an concentration impulse and the Fourier 
transforms of corresponding linear response func- 
tions. We observe how the input signal is transmit- 
ted through the pathway, with the effective pulse 
durations becoming progressively longer along the 
pathway (Fig. [3J left), including total receptor ac- 
tivity in a cell (A c ), phosphorylated kinase CheA, 
phosphorylated response regulator CheY, and fi- 
nally the motor, characterised by its probability of 
tumbling (Pew)- In Fig. [3J (middle) we show the 
corresponding linear response functions. 

The receptor activity acts as a high-pass filter, 
i.e. it transmits high-frequency signals, but not 
low-frequency signals. As can be seen from our 
simple model (cf. Eq. ([5])), this property is due to 
adaptation, which introduces the time-derivative of 
the signal Ac(t) up to the characteristic frequency 
u>m, eliminating the response to slowly changing at- 
tractant concentrations. The activity of chemore- 
ceptors is the input to further levels in the path- 
way. The response of CheA-P is fast, and shows 
no qualitative difference to the response of recep- 
tors in the frequency range shown. In contrast, due 
to the fast but finite rates of phosphorylation and 
dephosphorylation, preventing the CheY-P concen- 
tration to respond to rapidly changing input sig- 
nals, the response at the level of CheY is reduced 
at high frequencies. Similarly, the motor introduces 
another high-frequency filter due to slow switching 
between its two states. This additional filter can be 
deduced from Eq. ©, where the motor response 
function takes the response of chemoreceptors as 
input, and additionally introduces a characteristic 
cut-off frequency wp cw due to slow motor switch- 
ing rates. Hence, t he chemotaxis pat hway acts as 



a band-pass filter (Block ct all 119821 ). which only 



transmits input signals within a selected frequency 
range, which is of the order of 1 to 10 s. This time 
scale corresponds to the average time between two 
tumbles, allowing sensing of concentration changes 
during periods of running. 

As shown in Fig. [3J middle the phase tends to- 
wards 7r/2, i.e. a quarter peri od, at low frequencies . 



This has been analysed by Shimizu et al. ( 2010h 



only for the receptor complex activity. This phase 
difference is due to adaptation and represents the 
fact that the system takes the time derivative of the 
stimulus below the characteristic frequency ujm of 
adaptation. The phase shift of the receptor activity 
increases to 7r at high frequencies, indicating that 
the activity simply follows the output (a negative 
sign is due to the negative r esponse of the a ctivit y 
to attractant concentration; Shimizu et all 20101) . 
The phase at high frequencies for the response of 
CheA follows the phase of the receptor activity, ex- 
cept for a small increase of the phase shift. In con- 
trast, the phase of CheY and the motor increase 
significantly beyond 7r indicating that slow rates of 
modification and motor switching introduce a lag 
of the response behind the stimulus. 



Noise propagation 

To understand the noise characteristics of the mo- 
tor, we consider the noise sources and their trans- 
mission in the pathway. Each step in the signalling 
pathway is essentially probabilistic, hence, noisy: 
ligand diffusion and binding, receptor switching be- 
tween its functional on and off states, as well as re- 
ceptor methylation and demethylation, phosphory- 
lation and dephosphorylation of signalling proteins 
CheA, CheY and CheB, and switching of the rotary 
motor between its two states, CW and CCW rota- 
tion. To characterise fluctuations of the phosphory- 
lated signalling protein 8R(t) around its mean value 
(R(t)), we use the power spectrum Sr(u) and the 
variance (SR 2 ) = (R 2 (t)) - (R(t)) 2 (cf. Box 2). 

Analytical results for noise spectra. Consid- 
ering Eq. ([l])-((4|) with noise, we can analytically 
calculate power spectra (see Materials and Meth- 
ods). The power spectrum of activity fluctuations 
is given by 



S A » = N c - 



M 



+ LU 2 



(8) 

In this equation we considered fluctuations from 
receptor switching (first term in numerator), lig- 
and diffusion (second term), as well as the recep- 
tor methylation dynamics (third term) at each of 
the Nc receptor complexes per cell. We have as- 
sumed that fluctuations at different receptor com- 
plexes are independent. Therefore, we obtain the 
sum of Nc identical spectra for all complexes. The 
individual terms S a (u), S c (uj) and Qm are given by 
Eq. (HH, (TTTJ) and (JT^Jl in Materials and Methods. 
The frequency dependence of the ligand noise, as 
well as noise from receptor complex switching, in- 
dicates filtering of slowly varying fluctuations with 
frequencies below the characteristic frequency ujm 
due to adaptation. In contrast, only high-frequency 
fluctuations from the receptor methylation dynam- 
ics are filtered by the adaptation dynamics. This 
is due to finite rates of methylation and demethy- 
lation fluctuations introducing correlations in the 
receptor methylation level. 

The power spectrum of fluctuations in motor bias 
is obtained as 



- 0J 



(9) 



The first term represents transmitted noise from re- 
ceptor complexes, including the noise power spec- 
trum of the receptor activity and the sensitivity 
of motor switching rates to changes in activity uj 2 . 
The second term is motor switching noise. Both 
noises are filtered by the motor, as its finite rates 
of switching introduce correlations with character- 
istic frequency wp cw . The noise spectra of the full 
pathway including the phosphorylation reactions 
are shown in the SI. 
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Figure 3: Chemotaxis pathway response and noise spectra. (Left and Middle) Response upon impulse 
stimulation with attractant MeAsp. (Left) Time courses for MeAsp concentration c, total activity of 
receptors A c , CheA-P and CheY-P, and motor bias Pew- (Middle) Fourier transformed ligand signal, as 
well as response functions. Symbols correspond to numerical simulations and solid lines to analytically 
calculated response functions. (Right) Noise spectra of ligand and for the total activity of receptors A c , 
CheA-P and CheY-P, and motor bias Pew (thick black lines). Also shown are the contributions to 
the spectrum from ligand binding (L; thin solid black lines), receptor switching (RS; thin red lines 
and symbols), receptor methylation and demethylation (M; green lines), as well as phosphorylation and 
dephosphorylation of CheA, CheY and CheB (P; blue lines) and motor switching (MS; dashed grey line). 



Noise filtering along the pathway. In Fig. [3] 
(right), we show the power spectrum of fluctua- 
tions at the various levels of the signalling path- 
way, i.e. total receptor activity, CheA-P, CheY-P 
and the motor. We also plot the individual con- 
tributions from processes generating noise, namely 
ligand diffusion, receptor switching, methylation 
and demethylation of receptors, and phosphoryla- 
tion and dephosphorylation of proteins, as well as 
motor switching. This allows us to follow how noise 
is generated and transmitted at the various levels 
of the pathway. The noise spectrum of the recep- 
tor activity has its largest contribution at low fre- 
quencies, which originates in the receptor methy- 
lation and phosphorylation dynamics. Most of the 
fluctuations from phosphorylation stem from CheB 
(the separate contributions to the phosphorylation 
noise are not shown in Fig. [3J right). At high- 
frequencies, the activity noise spectrum is flat. This 
is due to ligand and receptor switching noise, which 
is removed at low frequencies by adaptation, but 
not at high-frequencies. The general behaviour of 
the noise spectrum corresponds to the simplified 
model (cf. Eq. ©). 

The noise spectrum of CheA-P has generally 
the same shape as the activity spectrum with a 



large low-frequency component, mainly due to re- 
ceptor methylation and CheB phosphorylation dy- 
namics. This spectrum also has an almost flat 
high-frequency behaviour in the frequency range 
shown. Apart from ligand and receptor switching 
noise, the flat part of the spectrum is largely deter- 
mined by fluctuations from CheA autophosphoryla- 
tion, which has roughly the same shape as activity 
noise at high frequencies because autophosphoryla- 
tion depends on the receptor activity. 

The noise spectrum of CheY-P is also largest at 
low frequencies. However, at high frequencies the 
spectrum falls off as noise is filtered due to the finite 
rates of CheY phosphorylation and dephosphoryla- 
tion, which introduce correlations in the fluctua- 
tions. 

The motor introduces another layer of filtering 
of transmitted noise with the characteristic motor 
switching frequency wp cw (cf. Eq. ©). Hence, 
transmitted noise is reduced by two filters in the 
frequency range shown, namely due to the CheY-P 
and motor dynamics. However, the main contribu- 
tion to the spectrum is due to the motor switching 
itself, which is reduced only by a first-order filter 
with characteristic frequency wp cw . 
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Cell-to-cell variation of motor 
behaviour 

How are the signal response, fluctuations and the 
signal-to-noise ratio (SNR) affected by changing 
parameters of the pathway such as size of recep- 
tor complexes, protein concentrations and reaction 
rate constants? In this section, we discuss the ef- 
fect of cell-to-cell variation on the power spectrum 
of the motor. In the next section, we discuss the 
SNR and its contributions, and how they depend 
on receptor complex size and adaptation rates. 

According to our model parameters obtained 
from fitting the Fourier transformed linear response 
to data, the main contribution to the power spec- 
trum comes from the steady-state switching of the 
motor between CCW and CW state. However, 
cell-to-cell variation in protein content and motor 
switching rates can lead to modifications of the 
largely Lorentzian-shapcd spectrum. These mod- 
ifications are caused by the transmitted noise from 
receptor methylation and phosphorylation dynam- 
ics (green and blue lines in Fig. [3] right). Specifi- 
cally, Fig. shows the motor power spectrum for 
increased motor switching rates as well as reduced 
adaptation rates and number of chemoreceptors in 
a cell. In all cases the low-frequency component of 
the transmitted noise becomes more prominent. 

An increased low-frequency component has been 
observed in the mo tor power spectrum for c ells with 
low motor bias by iKorobkova et al. | (|2004) . These 
authors measured the motor power spectrum for 
wild-type cells, as well as mutants lacking the sig- 
nalling pathway. The mutant's spectrum represents 
the component to the power spectrum from steady- 
state motor switching only. Wild-type cells showed 
a large low-frequency component compared to the 
mutants. Figure 0)3 shows that our model can re- 
produce these experimental data (shown in the In- 
set). 

Signal-to-noise ratio at the mo- 
tor 

To characterise how signals are transmitted in the 
presence of noise, we define the SNR at the level of 
the motor as 

A d2 

CW 



SNR = 



<«3w> 



(10) 



with AP3 W 



and {5 Pew ) defined in Box 2. For 
optimal signalling this ratio should be maximised. 
For simplicity, we only discuss the receptor activity 
in the text, while in the figures we additionally show 
the contribution from phosphorylation processes as 
transmitted to the motor. 

Optimal receptor complex size 

Receptor complexes amplify small signals propor- 
tionally to their size N. However, also concen- 
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Figure 4: Effect of cell-to-cell variation on power 
spectrum of the motor. (A) Predictions about how 
different cell parameters affect the motor power 
spectrum, in particular its low-frequency compo- 
nent. The black line (wild-type WT1) is the same 
as the total motor spectrum in Fig. [3] (right). The 
motor spectra for increased motor switching rates 
(MS; red line), as well as reduced rates of receptor 
methylation and demethylation (M; green line) and 
the total number of receptors (blue line) are shown 
as well. Arrows indicate the features in the spectra 
that are affected. Dashed box is the area shown 
in panel (B). (B) Motor spectrum of cells with low 
motor bias (black line), as well as the component 
from steady-state motor switching only (grey line) . 
(Inset) Measured power spectra for wild-type cells 
(WT2) with low motor bias (black) and mutant 
lacking the signalling pathways (grey). Axes are 
the same as in the main panel. Error bars indicate 
the measurement uncerta inty. Spectra wer e trace d 



from data presented by IKorobkova et all (|2004l) 



Model parameters are listed in the SI. 



tration fluctuations are expected to be amplified. 
Hence, we hypothesize that the receptor complex 
size could be optimised to yield a balance of ad- 
vantageous amplification of signals and detrimental 
amplification of input noise. 

In Fig. [5] (top) we show the integrated motor 
response APq W (see Box 2) to a step stimulus 
for varying background concentration and receptor 
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Figure 5: Varying ambient concentration and receptor complex size. ( Top) Integrated response of the 
motor bias. (Middle) Variance of the motor bias including only contributions from receptor switching, 
ligand diffusion, methylation and phosphorylation. The individual contributions are shown to the right 
of the main panel. (Bottom, left) SNR based on the signal response and variance shown in the top and 
middle panel. (Bottom, right) SNR as a function of receptor complex size at ambient concentration 0.02 
(solid), 0.03 (dashed) and 0.05 mM (dotted line). Optimal complex size is indicated by arrows. 



10 



complex size. We assume that the step stimulus 
size is a constant fraction of 10 percent of the back- 
ground concentration. The integrated response has 
a characteristic variation with background concen- 
tration with the maximum in the sensitivity range 
of Tar receptors (indicated by their dissociation 
constants). Furthermore, the response increases 
with receptor complex size N. We calculated the 
integrated signal response of the receptor activity 
(see Materials and Methods). This quantity scales 
linearly with receptor complex size, A^4^ oc N, due 
to coherent addition of the signalling responses of 
different receptor complexes, amplification of con- 
centration changes by receptor complexes, as well 
as filtering by adaptation. 

In the middle panel of Fig. [5j we show the vari- 
ance (i.e., the integrated noise power spectrum, see 
Box 2) of the transmitted noise of the pathway at 
the level of the motor. Only the contribution to 
the variance from ligand diffusion depends on the 
background concentration. Compared to the signal 
response, the maximum of the variance is shifted to 
a slightly lower concentration. The contribution to 
the variance from switching of receptor complexes 
is relatively small compared to the other contri- 
butions and roughly constant with receptor com- 
plex size, whereas those from ligand diffusion, re- 
ceptor methylation and phosphorylation dynamics 
increase with receptor complex size. 

To understand these behaviours of the variance 
more intuitively, we analysed the receptor activ- 
ity analytically (for details of the calculation, see 
Materials and Methods). We find, the contribu- 
tion to the variance of the receptor activity from 
receptor switching is indeed constant, independent 
of N. The contribution from ligand diffusion scales 
steeply as A^ 2 , the difference between ligand noise 
and ligand signal amplification being due to (i) 
noise from different complexes is added up incoher- 
ently, and (ii) the main contribution to the variance 
coming from high-frequency ligand noise, which is 
not filtered by adaptation. The contribution from 
receptor methylation grows approximately linearly 
with receptor complex size as a result of the inco- 
herent addition of fluctuations at different receptor 
complexes and the sensitivity of the receptor com- 
plex activity with respect to changes in methyla- 
tion level increasing proportionally with AT. The 
contribution to the variance from phosphorylation 
processes grows with receptor complex size simi- 
lar to the contribution from the methylation dy- 
namics. Overall, the total variance of transmitted 
noise at the level of the motor has contributions 
from receptor switching, the dynamics of receptor 
methylation, and phosphorylation. The latter is 
approximately constant or grows slower than the 
amplified signal response, whereas the component 
from ligand diffusion increases steeper than the sig- 
nal response with growing receptor complex size. 

The resulting SNR, i.e. the ratio of integrated 
signal response and variance of the noise, is shown 
in Fig. [5] (bottom). The SNR is largest at back- 



ground concentrations in the sensitivity range of 
the Tar receptor. Furthermore, due to the differ- 
ent dependencies of the signal and the noise on the 
receptor complex size, the SNR has a maximum 
at a particular receptor complex size (Fig. [5] bot- 
tom, right). The SNR grows below that complex 
size due to signal amplification, while the amplified 
ligand noise from ligand diffusion is still below the 
internal noise level from receptor switching and re- 
ceptor methylation and phosphorylation dynamics. 
Above the optimal receptor complex size, the SNR 
decreases because the ligand noise is amplified more 
than the signal. 

Optimal adaptation rates 

Adaptation filters slow input signals, its speed de- 
termining what input frequencies are transmitted 
by the pathway. Furthermore, the adaptation dy- 
namics filters input noise. Hence, adaptation rates 
may be expected to be optimised for signal and 
noise propagation. 

Figure [6] shows the integrated signal response at 
the level of the motor for varying rates of receptor 
methylation (7^) and dcmcthylation (7s). Vary- 
ing these parameters describes changing the con- 
centrations of receptor modification enzymes CheR 
and CheB. Interestingly, varying the two parame- 
ters has different effects on the signalling response: 
the integrated signal response increases for van- 
ishing 7^, whereas it decreases for vanishing 73. 
There are two effects that contribute to this be- 
haviour: Firstly, if the concentration of one of the 
receptor modification enzymes is reduced, the re- 
ceptors becomes modified predominantly by the op- 
posing enzyme, hence driving the receptor activ- 
ity towards saturation (A* = or A* = 1). This 
effect would tend to quench the response by re- 
ceptors. Secondly, as the enzyme concentration is 
reduced, adaptation times increase. Hence, this 
effect increases the integrated signal response as 
the time the receptor activity deviates from the 
adapted state increases. According to calculations 
shown in Materials and Methods for the integrated 
response of receptors, the first effect dominates in 
the case of reduced jb- Due to the strong activ- 
ity dependence of the demethylation rate, reducing 
the demethylation rate constant effects the adapted 
activity of receptors strongly. Hence, receptors are 
quickly driven into saturation for vanishing 75 . In 
contrast, in the case of reduced the second ef- 
fect dominates and the increased adaptation time 
leads to an increased integrated signal response. At 
large methylation and demethylation rates, adap- 
tation times are reduced leading to a decreasing 
integrated signal response. 

The variance of fluctuations is shown in the mid- 
dle panel of Fig. [SJ The individual contributions 
from transmitted noise at the level of the motor 
look qualitatively similar. All contributions de- 
crease both for vanishing 7jj and 75 consistent with 
calculations for the variance of the receptor activity 



11 



Signal [s] 




Methylation rate constant y 



Figure 6: Varying receptor methylation and demethylation rate constants 7^ and 7b, respectively. (Top, 
left) Integrated response of the motor bias. [Top, right) Illustration of the effects of vanishing 7^ and 
7b on adapted activity (indicated by dot and circle along dose-response curve; left), as well as on time 
courses (right) for three cases, 7_r — > (top right), jr/jb — const (middle right) and 7s — > (bottom 
right). For further explanation see text. (Middle, left) Variance of the motor bias including only from 
receptor switching, ligand diffusion, methylation and phosphorylation. The four individual contributions 
are shown to the right of the main panel. (Bottom, left) SNR based on the signal response and variance 
shown in the top and middle panel. 
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in Materials and Methods. 

The SNR is shown in the bottom panel of Fig. [5J 
The SNR increases for vanishing 7^ and decreases 
for vanishing 73. According to Fig.EJ a large SNR 
is obtained for small and large jb, correspond- 
ing to the parameters of our model. 

Fluctuation-response relation- 
ships 



0.005 r 



IPark et all fl2010h presented the idea that the sig- 
nalling response to concentration signals and fluc- 
tuations in the chemotaxis pathway are not inde- 
pendent of each other, because they are produced 
by the same molecular interactions. Specifically, 
based on measurements at the level of the motor 
these authors proposed a fluctuation-response the- 
orem, namely an approximate linear relationship 
between the adaptation time to step stimuli and 
the variance of fluctuations in CheY-P concentra- 
tion. 

Using our model, we tested this hypothesis and 
varied the adaptation rates, as well as the to- 
tal CheY concentration in a cell, resulting in a 
shifted adapted CheY-P concentration at steady 
state. We find that the variance of CheY-P (nor- 
malised by the squared adapted value) decreases as 
the adapted CheY-P value increases except for very 
small adapted CheY-P concentrations (Fig. [JjA.), 
indicating that the relative strength of fluctua- 
tions decreases as expected. In Fig. [7)3 we show 
the adaptation time, approximated by the inverse 
of the characteristic frequency due to adaptation, 
plotted against the variance of CheY-P. We find 
that at low adaptation times (thick line styles of 
the curves) , the adaptation time increases with the 
variance of CheY-P, indicating that cells with large 
fluctuations also respond longer to concentration 
signals. In contrast at long adaptation times, the 
adaptation time decreases with increasing variance 
of the pathway (grey parts of the curves). This 
behaviour can be directly traced back to the non- 
monotonic variance shown in Fig. [BJ It is maximal 
when the adapted CheY-P concentration is about 5 
/xM, i.e. when typically half of CheY is phosphory- 
lated. The exact relationship depends on what pa- 
rameteris varied, exemplified by the different curves 
in Fig. [7)3. For each parameter and small adapta- 
tion times, we fi nd an approximat e linear relation- 
ship in line with Park et aTl f|2010h . see Inset. 



Discussion 

Biological signalling pathways employ biochemi- 
cal reaction networks and molecular state transi- 
tions to sense and process signals from the environ- 
ment. Fluctuations inherent in these processes de- 
termine the signals which can reliably be transmit- 
ted. Here, we studied the signal and noise propaga- 
tion in the E. coli chemotaxis signalling pathway, 



0.004 



0.003 - 



Q_ 0.002 

> 

o 0.001 



100 



Varied parameter: 

— r„ 
-- y B 
■■■■ m„ 



4 6 

[CheY-P]* [nM] 



10 



80 



CD 

■I 60 

c 

o 

I 40 
03 
"D 
< 



20 



Park et al. (2010) 




0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 

CheY-P noise [|aM 2 ] 



Figure 7: Fluctuation-response relationships. (A) 
Variance of CheY-P fluctuations (normalised by the 
squared adapted CheY-P value) as a function of 
adapted CheY-P value for varying adaptation rates 
7_r (solid line) and jb (dashed line) and total CheY 
concentration (dotted line). (B) Adaptation time, 
calculated as (XidA/dM) -1 , plotted against the 
variance of CheY-P fluctuations. Adapted CheY-P 
concentration of 5 fiM is indicated by the circles. 
Legend is the same as in (A). Our parameters are 
given in the SI. (Inset) Variance of CheY-P con- 
centrations plotted against the adaptation time as 
extracted f r om m easurements of motor rotation by 
Park et all (|2010h . The line is a linear fit through 
the data. 



which controls the bacterial swimming behaviour 
in chemical gradients. Specifically, we considered 
the dynamics of ligand diffusion, receptor methyla- 
tion and demethylation, receptor complex switch- 
ing between on and off, phosphorylation and de- 
phosphorylation of the kinase CheA, and response 
regulators CheY and CheB, as well as from ro- 
tary motor switching between CW and CCW di- 
rection. We assume cooperative chemoreceptor sig- 
nalling complexes, whose activities depend on lig- 
and concentrations and receptor methylation level , 
descr i bed by the MWC model (ISouriik and Bera. 



Endres and Wingreenl 20061 : lEndres et al 



2004 iMello and Tul. liooj iKevmer et al 



Clausznitzer et al. . 201dh 



2006; 



2008; 



We formulated a model which includes all pro- 
cesses in the signalling pathway, discussed in the 
SI. Not included is the dynamics of gene expression, 
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which is assumed to be much slower than the dy- 
namics of signalling processes. To make results in- 
tuitive we also presented a simplified version of the 
model, which only includes the dynamics of the ac- 
tivity of chemoreceptors, ligand concentration and 
receptor methylation level, as well as the motor dy- 
namics. To calibrate the model, we first collected 
experimental data sets for the signa l ling pathway 



and rotary motor (Block et all Il982t iSeeall et al. 
19861: IShimizu et all bOlolT and the motor switch- 



ing behaviour ( Korobkova et al. . 2004 . 2006). Us 
ing the Fourier transformed linear response func- 
tion, we subsequently fitted our model parameters. 
We found a range of parameters fitting different 
data sets, revealing a striking experimental varia- 
tion, which may require further characterization in 
the future. 

Despite the fitting, there is a discrepancy of our 
respo nse function and th e dat a at large frequen 



cies. block et al.l f|l982h and ISeeall et all (|l986h 



conjectured that the pathway is a third-order low- 
pass filter. In contrast, we find that the only rel- 
evant filters in that frequency range are due to 
CheY-P and motor dynamics, leading to only a 
second-order filter. One explanation for the miss- 
ing filter is that experimental concentration pulses 
were not short enough, leaving a signature from 
the input signal at large frequencies. Alterna- 
tively, additional processes such as a slow release of 
CheY-P fr om the chemosen sory complexes as dis- 
cussed by Blat et all (|l998l ) could lead to an ad- 
ditional filter. However, CheY-P/CheZ complex 
formation and pote ntial oligomerisation of CheY - 
P/CheZ comple xes ( Blat and Eisenbachl Il996bllal ; 
Eisenbachl 120041) are not expected to contribute to 
high-frequency filtering (SI). 

The motor behaviour is the final cell output, 
which contains characteristic noise signatures of all 
upstream signalling components, including the re- 
ceptors. We found that motor switching is the 
dominant contribution to the spectrum of the fluc- 
tuations in motor bias. This may be not surpris- 
ing as motor switching enables E. coli to tumble 
and change its swimming direction, and is there- 
fore crucial for its search strategy. However, low- 
frequency contributions from the dynamics of re- 
ceptor methylation and phosphorylation processes 
may be dominant in particular cells (Fig. [4]). We 
predict that due to cell-to-cell variation of pro- 
tein contents or fast motor switching, these low- 
freque ncy components b e come important. Specif- 



ically, iKorobkova et al ( 20061 ) measured power 
spectra in cells with low motor bias and found an 
increased low-frequency component as compared to 
mutant lacking the signalling pathway. Our model 
is able to reproduce these spectra, long correla- 
tions in motor bias may enable subpopulations of 
cells to swim further without tumbling or to tumble 
more frequently. 

Although chemotaxis is one of many capabilities 
a cell has and may not be optimised in isolation 
without the rest of the cell, we speculate the cell 



aims to maximise the SNR for most efficient sig- 
nalling and chemotaxis. We found that the SNR 
is maximised at particular receptor complex sizes 
similar to values of receptor cooperati yity extracted 



from FRET dose-response curves (jEndres et al 
2008). In line with the data, the "optimal" complex 
size increases with external ligand concentration, 
and hence with receptor methylation level. While 
our complex sizes appear overestimated, noise from 
ligand m olecules rebinding to th e same receptor 



complex (Bial ek and Setaveshgarl 120051 ) has not 
been considered here. This may well increase the 
noise level from external sources and hence decrease 
the predicted optimal receptor complex size. Using 
our model, we also analysed the effect of varying 
the methylation and demethylation rate constants. 
We found that a large SNR is obtained for small 
methylation and large demethylation rate constant, 
corresponding to our fitted model p arameters from 
FRE T dose-response curves (jClausznitzer et al, 
20101) . 



We predict that the noise power spectrum of the 
motor reveals a significant low-frequency compo- 
nent when transmitted pathway noise becomes im- 
portant compared to motor switching noise (Fig. 0]) 
To test this prediction, adaptation rates can be 
varied using cells expressing different amounts of 
CheR and CheB from an inducible plasmid. Al- 
ternatively, the natural variability in protein ex- 
pression between cells can be exploited. Numerous 
wild-type cells could be monitored. By extracting 
the adaptation times for chemotactic stimuli, the 
adaptation rate constants for individual cells can 
be inferred. Subsequently, the same cells can be 
used to measure long time traces of motor switch- 
ing and noise spectra can be calculated. 

Our full-pathway model allows us to connect to 
a large variety of data and literature. For instance, 
we also investigated how the response to concentra- 
tion signals is related to fluctuations in the chemo- 
ta xis pathway (Fig. [71 ), similar to ideas presented 
by IPark et all (|2010l) . We do not find one unify- 
ing fluctuation-response theorem, but a nonlinear 
trend for the relationship between adaptation times 
and variance in CheY-P. For small adaptation times 
we fi nd an approxima te linear relationship in line 
with lPark etaD (|2010h . 

To describe the noise, the langevin approxi- 
mation is expected to work for the phosphoryla- 
tion and dephosphorylation of the abundant pro- 
tein CheY. However, its applicability is less clear 
for receptor signalling due to both extrinsic lig- 
and noise and intrinsic noise from receptor methy- 
lation. Furthermore, the switching of the binary 
motor may constitute relatively large noise. As 
shown in the SI, the Master equation and langevin 
approxima tion yield the same results for receptor 



signalling (jAauino et all 120111) . As for the motor 



we explicitly tested that the statistical properties 
of the time series obtained for two-state switching 
and langevin equation are the same. For constant 
rates, as well as for noisy rates due to fluctuations 
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Box 3 Comparison of E. coli chemotaxis and other two-component systems. 




In the chemotaxis pathway (left panel) and other two- 
component systems (right panel), a sensor kinase is acti- 
vated by a cell-external signal, upon which it autophos- 
phorylates and passes on a phosphoryl group to its re- 
sponse regulator, which typically induces a transcrip- 
tional response. The time course of the final output, 
i.e. gene expression, can be directly mapped onto the 
binary output of the chemotaxis pathway. 

There are particularly well-studied examples of two- 
component systems: (i) the VanS (kinase) /VanR (re- 
sponse regulator) system conferring vancomycin re- 
sistan ce in Gram positive bacteria dHutchings et all 
2006 ), (it) quorum sensing in Vibrio harveyi, where 
the three kinases LuxN, LuxQ and CqsS respond 



to different autoinducers and first phosphorylate the 
phosphotransferase LuxO (which has no equivalent 
in the chemotaxis pathway), whi ch then phosphory- 
lates the response regulator LuxU l|Henke and Basslerl . 
2004), and (lii) the phosphorelay controlling sporula- 
tion in Bacillus subtihs. The relay contains at least 
four kinases KinA-KinB and the phosphotransferase 
SpoOF, which phosphory lates the response regulator 
SpoOA l| Jiang et al.l . l2000h . 

Most other bacterial chemotaxis pathways are more 
complex than E. colts pathway. For instance, the pho- 
tosynthetic bacterium Rhodobacter sphaeroides has sev- 
eral homologues of each of the chemotaxis proteins in 
E. coli (|Porter et all l2008h . Interestingly, Rhodobac- 
ter has two chemotaxis receptor clusters, one polar 
cluster similar to E. coli and one cytoplasmic clus- 
ter, which is thought to sense the metabolic state of 
the cell. Both clusters need to be present for chemo- 
taxis. The soil bacterium Bacillus subtilis has three 
adaptation systems |Glckas et al., 2010j): one based on 
methylation and demethylation of receptors similar to 
E. coli and two independent of receptor methylation, 
the CheC/CheD system and the CheV system. Fur- 
thermore, in B. subtilis, sensory adaptation is not de- 
termined by the level of receptor methylation but the 
location of methylation at receptors. 



in CheY-P concentration, the power spectra ob- 
tained for the two processes are the same (SI). In 
the SI, we further show that noise terms are indeed 
sufficiently small that linearisation of the pathway 
equations is justified. 

The bacterial chemotaxis pathway is a mem- 
ber of the large class of two-component sys- 
tems, containing hundreds of closely related path- 



response, 



virulence and 



ways involved in stress 

inter-cell communication ( Laub and GoulianL 2007 
Souriik and Armitagel . l2010h . In these pathways, 
activation of a sensor histidine kinase results in 
its autophosphorylation, and subsequently in phos- 
phorylation of a response regulator, which typi- 
cally binds to DNA and regulates gene expres- 
sion (Box 3). The final output, i.e. activation of 
gene expression, is again binary and hence similar 
to the bacterial chemotaxis pathway. The analysis 
presented here may also help elucidate the design 
of many other pathways and clarify the computa- 
tional problems cells try to solve. 

Materials and Methods 

MWC model for activity of receptor 
complexes 

The MWC model describes signalling by receptor com- 
plexes, which can switch between their on and an off 
states. The average activity of a complex is given by 



A = 



N 



1 + e F (cM) ■' 



(11) 



ranging from zero to N. The free-energy difference 
F(c, M) between the on and off state is 



F(c,M) = N ■ 



-M + N 



In 



c/Kf 



1 + c/Kf 



+v s In 



1 + c/K° 



(12) 



l + c/Kf / 

which is a function of the concentration c present at the 
receptor complex site and the methylation level M of 
the receptor complex. The methylation level of a com- 
plex is the sum of methylation levels of all receptors 
in a complex. Here, we consider two receptor types, 
Tar (indicated by index a) with fraction v a of recep- 
tors in the complex, and Tsr (indicated by index s) 
with fraction v a of receptors. Receptors are sensitive 
to attractant MeAsp with dissociation constants K° n 
and K° B in the on and off state, respectively. We use 
the following parameters for the MWC model for re- 
ceptor complexes: K°f = 0.02 mM, K° u = 0.5 mM 
Kf = 100 mM and Jf ° n = 10 6 mM i|Kevmer et all 
120061 : IClausznitzer et all l20ld ) . 

Noise sources 

Switching noise. The switching noise T]A(t) in 
Eq. (ffj is due to the switching of each receptor complex 
between on and off states. We assume the switching to 
be a fast process, which can be described by the fol- 
lowing dynamics for the the probability of a receptor 
complex to be on, a: 

^=ki(N-a)-k 2 a + ria(t). (13) 
at 

The noise term rj a (i) is a Gaussian white noise with 
zero mean and noise intensity Q a = 2k2A* , where we 
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used that the receptor complex activity A — (a) which 
is equal to the (quasi) steady-state activity of a, and 
A — A* when adapted. The power spectrum of a due 
to switching between on and off states is 



Motor switching noise. The noise in motor 
switching rate in Eq. © is assumed to be a Gaus- 
sian white noise term with zero mean, noise intensity 
Qp cw = 2k+(A*)(l — Pqw) an d power spectrum 



S a (uj) = 



03 2 + (fcl + fc 2 ) 2 



(14) 



where ki + k 2 is the characteristic frequency of switch- 
ing. Hence, the high-frequency component of fluctu- 
ations 5a(t) is reduced due to averaging by the finite 
rates of switching. Hence, the power spectrum of ac- 
tivity fluctuations r)A(t) is 



(15) 



Ligand noise. The number of ligand molecules in 
the vicinity of a receptor complex fluctuates due to 
binding/unbinding, and potential rebinding of pre- 
viously bound ligand molecules at this complex , 
as well as diffusion jBialek and Setaveshgarl . 120051 : 
lEndres and Wingreenl. 120091 1. Here, we use a simpli- 
fied description of diffusion to calculate the spectrum of 
noise in the ligand dynamics rj c (t) in Eq. Consider 
a volume whose dimensions are given by the diameter 
of a receptor complex s — V Nsr, where s r = 1 nm is 
the size of a receptor dimer (|Hazelbauerl.ll992h . The 
change of ligand-molecule number L in this volume is 
determin ed by the ex change rate ko ~ D/(2s 2 ) due to 
diffusion (|Berd . ll993r i: 



dL . 3 



L)+r !L (t) 



(16) 



where koL is the rate of molecules moving out of the 
volume by diffusion, and ko times the mean concentra- 
tion Co in solution serves as a proxy of the rate of lig- 
and molecules moving into the volume. The noise term 
r)L(t) is assumed to be Gaussian and white, with zero 
mean and noise intensity Ql = Dsco. The power spec- 
trum of the number L and concentration c of molecules 
at receptor complex j is 



S L (u) 



Dscq 



S c (u) 



S L (u) 



(17) 



where s 6 is the squared volume given by the dimen- 
sion of the receptor complex. The zero-frequency limit 
of the power spectrum of the ligand concentration 
S c (0) = co/(Ds), which corresponds to cal culations by 
Berg and Purcell (iBerg and Purceli . 1977) and Bialek 
and Setayeshgar (|Bialek and Setaveshgar , 2005I - ) for the 
uncertainty in sensing ligand concentration. The noise 
rye (i) in Eq. @ is related to rate of change of the lig- 
and concentration, similar to the considerations of the 
switching noise above. Hence, the power spectrum of 
the ligand fluctuations rj c (t) is 



S Vc (uj) = UJ 2 Sc{u). 



(18) 



Methylation noise. The size of fluctuations in the 
rate of methylation of a receptor complex j in Eq. © 
is estimated from the average rates of methylation and 
demethylation at the adapted state, respectively. The 
noise TjM(t) is assumed to be Gaussian and white, with 
zero mean, noise intensity Qm ~ 2'Jr(N — A*) and 
power spectrum 



SvaA^) = Qm- 



Calculation of response functions 



(20) 



After linearising around the steady state and insert- 
ing the Fourier transforms we obtain for the simplified 
model 



- iujAA = 

—iuiAM - 
-icuAPcw — 



■ dA . - 
~ tUJ 'dM 



. dA . „ 
tuj-^—Ac 
oc 



(21) 

-cjiAA (22) 
w 2 Ai c -a;p cw APcw, (23) 



where 



uii = 7R + 3 7s A* 2 = 7/i (3 - 2A* r )/A* r 



(24) 



with A* = N ■ A* r « N/3 (jSouriik and Berd . l2002h 
the adapted activity of a receptor complex, A* denot- 
ing the adapted activity of individual receptors. In 
the second equality we have used that at the adapted 
state 7_r(-/V — A*) = jbA* 3 . The parameter 0J2 = 



(1 — Pqw) ~qa~ ~ Pcw~g~fi~~ ls t ne derivative of the motor 
switching rates with respect to activity, and wp cw = 
+ is a characteristic frequency due to motor 
switching at steady state. AA is the response of ev- 
ery receptor signalling complex, and AA C = NqAA is 
the activity response of all receptor complexes in a cell. 
Solving for AA C and APcw, and division by the stim- 
ulus Ac yields the response functions in Eq. ([5j) and 

Calculation of noise power spectra 

To calculate spectra, we linearise the deterministic 
parts of Eq. (UJ-© similar to the calculation of the 
response functions, and formally Fourier transform the 
equations. We obtain 



■jSA c 



jSMi 



(25) 



-uiSAj + fj M:l (26) 
uj 2 5A c - cjp cw <5Pcw + fjPcw (27) 



We solve for the Fourier transformed activity fluctua- 
tions 8 A c and obtain 



8A C 



mh + + Ej va 3 



(28) 



which yields the power spectrum in Eq. l[5). The pa- 
rameter ujm ~ uj\dA/dM, and we used Eq. (|15f) and 

From Eq. (|27|l we obtain for the Fourier transformed 
fluctuations in the probability of tumbling mode SPcw 



SPc^ 



uj 2 SA c + t)p cw 



LO Pa 



(29) 



(19) and their power spectrum is given by Eq. 
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Integrated signal response, variance 
and SNR 

Optimal receptor complex size. The inte- 
grated response of the receptor activity to a step stim- 
ulus is 



AAi 



dw|XA c (w)Ac(cj)|^ 



(30) 



where we inserted Eq. (T5J. Hence, the activity response 
scales as AA 2 oc (N tot /N) 2 (TV 2 ) 2 /TV oc TV, where we 
used that Nc = Ntot/N with TV to t the total number of 
receptors in a cell. 

The variance of the receptor activity is given by the 
integral over the power spectrum of activity fluctua- 
tions Eq. iJH 



(SA 2 C ) = 



Nc r' 1 ^[Sa{u) + (%) 2 S c { U j\ 
— — / duj = - 



2tt 

+ 



CO 2 + UJ 2 M 



Nc r 

27T J_ T -X 



c 2 + o, 2 ■ 



where we consider the frequency range relevant for mo- 
tor switching indicated by r _1 ~ 0.1 ... 1 Hz. 
The contribution from receptor switching is 



(6Al) a 



Nc_ 
2n 



UJ 2 S a (uj) 



2k 2 AtN u 



+ u)\ { 7rr(fci + fc 2 ) 2 



(32) 

where we used Q a and inserted Eq. (|14[1 for the power 
spectrum of receptor switching noise and used that it 
is almost constant and equal to its zero-frequency value 
over the integration range. Furthermore, the factor 
w 2 /(oj 2 + ujjj) ss 1 and A* — A* /N is the adapted 
activity of an individual receptor. Hence, according to 
this simple calculation the contribution to the variance 
from receptor switching is roughly constant with recep- 
tor complex size. 

The contribution from ligand diffusion is 



(SA 2 C ) C = 



TVc 
2tt 



(lid- 



LU 2 +Uji 



(33) 



{SA 2 )m grows approximately linearly with receptor 
complex size. 

The SNR grows linearly with TV for small complex 
sizes, and decreases as TV -2 for larger complex sizes, 
resulting in an optimal medium receptor complex size, 
in qualitative agreement with Fig. [5] 

Optimal adaptation rates. The integrated sig- 
nal response of the receptor activity Eq. (|30[) 



AA 2 C = 



8A 
dU 



(35) 



where the numerator expresses the initial response of 
receptors of concentration changes and the denom- 
inator the filtering by adaptation. The sensitivity 
dA/dc = NA*(1 - A*)h(c), where h(c) = dF/dc, 
wi = 7b+37s/V 2 (A;) 2 , and dA/dM = NA*(l-A*)/2. 
The adapted activity can be obtained analytically for 
our simplified model from the steady state of the methy- 
lation dynamics Eq. ©, 



(31) A r * 



3 1 



4 27 



? i, 



(36) 

and is only a function of the ratio /3 = Jr/jb- Expand- 
ing the adapted activity around A* = (for 7_r — > 0) 
yields A* oc 7^/ 3 , and around A* = 1 (for jb — > 0) 
yields A* oc 7_g. Similarly, ui\ oc r y 2 ^ 3 (const. +7^/ 3 ). 
Hence, dA/dc oc 7^ 3 (7s) and u)\dA/dM oc ■y'^ 3 (73). 

The initial response to concentration changes de- 
creases slower than adaptation times, resulting in an 
increased signal response for vanishing jr. For vanish- 
ing 7s, the initial response to concentration changes 
decreases faster than adaptation speed, hence yielding 
a vanishing signal response. The overall dependence of 
the integrated signal response is A A 2 oc 7^ ly ^ (7s) for 
7_r — > 0(73 — > 0). For the contributions to the vari- 
ance of the receptor activity from receptor switching, 
ligand diffusion and receptor methylation dynamics we 



obtain {8A 2 ) a oc j L R /3 ( 7s ), {SA 2 ) C oc 7^ ( 7 |) and 
(SA 2 ) M oc j R (7 2 / 3 ), respectively. 

Hence, according to our simplified model the SNR of 
the receptor activity goes as SNR oc J R 2 ^ 3 (7f/ 3 ), in 
qualitative agreement with Fig. [5] 



where {8c 2 } — co/(Dst) is the variance of the ligand 
concentration measured during the time interval r. We 
used Eq. ()17[) and the same argument as for the switch- 
ing noise to calculate the integral. Hence, the contribu- 
tion to the variance from the ligand diffusion grows as 
{8A 2 ) C oc TV 3 as a result of incoherent addition of noise 
from different receptor complexes and the sensitivity 
dA/dc increasing as TV 2 . 

The contribution to the variance from receptor 
methylation is 



2* {dMj J_ T -r ^ + 



Qm 



2TV tot 7 fl (l - A* r )A* r dA 
ui 8M 



Acknowledgements 

We thank Richard Berry, Martin Buck, Tom Duke 
and William Ryu for helpful discussions. RGE was 
supported by Biotechnological and Biological Sci- 
ences Research Council grant BB/G000131/1 and 
the Centre for Integrative Systems Biology at Im- 
perial College (CISBIC). 



Author contributions 



where we defined ui = Jr + 37s TV 2 (A*) 2 , inserted 
Qm = 27 fl TV(l - A*) and u M - •'■>''••" 



(34) DC and RGE conceived and designed the study, 
performed analytical calculations, analysed the 
data, and wrote the paper. DC performed com- 
= uj\(dA/dM). Hence, puter simulations. 



17 



References 

Acar, M., Pando, B. F., Arnold, F. H., Elowitz, 
M. B. and van Oudenaardcn, A. (2010). A Gen- 
eral Mechanism for Network-Dosage Compensa- 
tion in Gene Circuits., Science 329: 1656-1660. 

Aquino, G., Clausznitzer, D., Tollis, S. and Endres, 
R. G. (2011). Optimal receptor-cluster size de- 
termined by intrinsic and extrinsic noise., Phys 
Rev E 83: 021914. 

Baker, M. D., Wolanin, P. M. and Stock, J. B. 
(2006). Systems biology of bacterial chemotaxis., 
Curr Opin Microbiol 9: 187-192. 

Berg, H. C. (1993). Random Walks in Biology, 
Princeton University Press, Princeton. 

Berg, H. C. (2000). Motile behavior of bacteria., 
Phys Today 53: 24-29. 

Berg, H. C. and Purcell, E. M. (1977). Physics of 
chemoreception, Biophys J 20: 193-219. 

Bialek, W. and Setayeshgar, S. (2005). Physical 
limits to biochemical signaling., Proc Natl Acad 
Sci USA 102: 10040-10045. 

Blat, Y. and Eisenbach, M. (1996a). Mu- 
tants with defective phosphatase activity show 
no phosphorylation-dependent oligomerization of 
CheZ, J Biol Chem 271: 1232-1236. 

Blat, Y. and Eisenbach, M. (1996b). Oligomeriza- 
tion of the phosphatase CheZ upon interaction 
with the phosphorylated form of CheY, J Biol 
Chem 271: 1226-1231. 

Blat, Y., Gillespie, B., Brcn, A., Dahlquist, F. W. 
and Eisenbach, M. (1998). Regulation of phos- 
phatase activity in bacterial chemotaxis, J Mol 
Biol 284: 1191-1199. 

Block, S. M., Segall, J. E. and Berg, H. C. (1982). 
Impulse responses in bacterial chemotaxis., Cell 
31: 215-226. 

Block, S. M., Segall, J. E. and Berg, H. C. (1983). 
Adaptation kinetics in bacterial chemotaxis., J 
Bactemol 154: 312-323. 

Clausznitzer, D., Oleksiuk, O., L0vdok, Sourjik, 
V. and Endres, R. G. (2010). Chemotactic re- 
sponse and adaptation dynamics in Escherichia 
coll., PLoS Comput Biol 6: el000784. 

Colman-Lerner, A., Gordon, A., Serra, E., Chin, 
T., Resnekov, O., Endy, D., Pesce, C. G. and 
Brent, R. (2005). Regulated cell-to-cell variation 
in a cell- fate decision system., Nature 437: 699- 
706. 

Detwilcr, P. B., Ramanathan, S., Scngupta, A. and 
Shraiman, B. I. (2000). Engineering aspects of 
enzymatic signal transduction: photoreceptors in 
the retina., Biophys J 79: 2801-2817. 



Eisenbach, M. (2004). Bacterial chemotaxis, in 
M. Eisenbach (cd.), Chemotaxis, Imperial Col- 
lege Press, pp. 53-215. 

Eldar, A. and Elowitz, M. B. (2010). Func- 
tional roles for noise in genetic circuits., Nature 
467: 167-173. 

Elowitz, M. B., Levine, A. J., Siggia, E. D. and 
Swain, P. S. (2002). Stochastic gene expression 
in a single cell., Science 297: 1183-1186. 

Emonet, T. and Cluzel, P. (2008). Relationship be- 
tween cellular response and behavioral variability 
in bacterial chemotaxis., Proc Natl Acad Sci U S 
A 105: 3304-3309. 

Endres, R. G., Oleksiuk, O., Hansen, C. H., Mcir, 
Y., Sourjik, V. and Wingreen, N. S. (2008). Vari- 
able sizes of Escherichia coli chemorcccptor sig- 
naling teams., Mol Syst Biol 4: 211. 

Endres, R. G. and Wingreen, N. S. (2006). Precise 
adaptation in bacterial chemotaxis through "as- 
sistance neighborhoods"., Proc Natl Acad Sci U 
S A 103: 13040-13044. 

Endres, R. G. and Wingreen, N. S. (2009). Ac- 
curacy of direct gradient sensing by cell-surface 
receptors, Prog Biophys Mol Biol 100: 33-39. 

Falke, J. J. and Hazelbauer, G. L. (2001). Trans- 
membrane signaling in bacterial chemorecep- 
tors., Trends Biochem Sci 26: 257-265. 

Glekas, G. D., Cates, J. R., Cohen, T. M., Rao, C. 
and Ordal, G. W. (2010). Site-specific methyla- 
tion in Bacillus subtilis chemotaxis: The effect of 
covalent modifications to the chemotaxis recep- 
tor McpB., Microbiology . 

Hazelbauer, G. L. (1992). Bacterial chemorecep- 
tors., Curr Opin Struct Biol 2: 505-510. 

Henke, J. M. and Basslcr, B. L. (2004). Three 
parallel quorum-sensing systems regulate gene 
expression in Vibrio harveyi., J Bacteriol 
186: 690214. 

Hornung, G. and Barkai, N. (2008). Noise propa- 
gation and signaling sensitivity in biological net- 
works: A role for positive feedback, PLoS Com- 
put Biol 4: e8. 

Hutchings, M. I., Hong, H. J. and Buttncr, M. J. 
(2006). The vancomycin resistance VanRS two- 
component signal transduction system of Strep- 
tomyces coelicolor., Mol Microbiol 59: 92335. 

Jiang, M., Shao, W., Perego, M. and Hoch, J. A. 
(2000). Multiple histidinc kinases regulate entry 
into stationary phase and sporulation in Bacillus 
subtilis., Mol Microbiol 38: 535-342. 

Kcymer, J. E., Endres, R. G., Skogc, M., Mcir, Y. 
and Wingreen, N. S. (2006). Chemosensing in 
Escherichia coli: two regimes of two-state recep- 
tors., Proc Natl Acad Sci U S A 103: 1786-1791. 



18 



Korobkova, E., Emonet, T., Park, H. and Cluzel, 
P. (2006). Hidden stochastic nature of a single 
bacterial motor., Phys Rev Letters 96: 58105. 

Korobkova, E., Emonet, T., Vilar, J., Shimizu, T. 
and Cluzel, P. (2004). From molecular noise to 
behavioral variability in a single bacterium., Na- 
ture 428: 574-578. 

Kubo, R. (1957). Statistical mechanical theory of 
irreversible processes i, Journal of the Physical 
Society of Japan 12: 570 - 586. 

Laub, M. T. and Goulian, M. (2007). Specificity 
in two-component signal transduction pathways., 
Ann Rev Genet 41: 121-145. 

Lestas, I., Vinnicombe, G. and Paulsson, J. (2010). 
Fundamental limits on the suppression of molec- 
ular fluctuations., Nature 467: 174-178. 

Meacci, G. and Tu, Y. (2009). Dynamics of the 
bacterial flagellar motor with multiple stators., 
Proc Natl Acad Sci U S A 106: 3746-3751. 

Mehta, P., Goyal, S. and Wingreen, N. S. (2008). 
A quantitative comparison of srna-based and 
protein-based gene regulation., Mol Syst Biol 
4: 221. 

Mello, B. A. and Tu, Y. (2005). An allosteric 
model for heterogeneous receptor complexes: un- 
derstanding bacterial chemotaxis responses to 
multiple stimuli., Proc Natl Acad Sci USA 
102: 17354-17359. 

Mora, T. and Wingreen, N. S. (2010). Limits of 
sensing temporal concentration changes by single 
cells., Phys Rev Lett 104: 248101. 

Mora, T., Yu, H., Sowa, Y. and Wingreen, N. S. 

(2009) . Steps in the bacterial flagellar motor., 
PLoS Comput Biol 5: el000540. 

Mora, T., Yu, H. and Wingreen, N. S. (2009). Mod- 
eling torque versus speed, shot noise, and rota- 
tional diffusion of the bacterial flagellar motor., 
Phys Rev Lett 103: 248102. 

Mortimer, D., Pujic, Z., Vaughan, T., Thompson, 
A. W., Feldner, J., Vetter, I. and Goodhill, G. J. 

(2010) . Axon guidance by growth-rate modula- 
tion., Proc Natl Acad Sci U S A 107: 52025207. 

Morton-Firth, C. J. and Bray, D. (1998). Predict- 
ing temporal fluctuations in an intracellular sig- 
nalling pathway, J Theor Biol 192: 117-128. 

Park, H., Pontius, W., Guet, C. C, Marko, J. F., 
Emonet, T. and Cluzel, P. (2010). Interdepen- 
dence of behavioural variability and response to 
small stimuli in bacteria., Nature 468: 819-823. 

Paulsson, J. (2005). Models of stochastic gene ex- 
pression, Phys Life Rev. 2: 157-175. 



Pedraza, J. M. and van Oudenaarden, A. (2005). 
Noise propagation in gene networks., Science 
307: 1965-1969. 

Porter, S. L., Wadhams, G. H. and Armitage, J. P. 
(2008). Rhodobacter sphaeroides: complexity in 
chemotactic signalling, Trends in Microbiology 
16: 251 - 260. 

Scharf, B. E., Fahrner, K. A., Turner, L. and Berg, 
H. C. (1998). Control of direction of flagellar 
rotation in bacterial chemotaxis., Proc Natl Acad 
Sci USA 95: 201-206. 

Segall, J. E., Block, S. M. and Berg, H. C. (1986). 
Temporal comparisons in bacterial chemotaxis., 
Proc Natl Acad Sci U S A 83: 8987-8991. 

Shimizu, T. S., Tu, Y. and Berg, H. C. (2010). A 
modular gradient-sensing network for chemotaxis 
in Escherichia coli revealed by responses to time- 
varying stimuli., Mol Syst Biol 6: 382. 

Sourjik, V. (2004). Receptor clustering and signal 
processing in E. coli chemotaxis., Trends Micro- 
biol 12: 569-576. 

Sourjik, V. and Armitage, J. P. (2010). Spatial 
organization in bacterial chemotaxis., EMBO J 
29: 2724-2733. 

Sourjik, V. and Berg, H. C. (2002). Receptor sen- 
sitivity in bacterial chemotaxis., Proc Natl Acad 
Sci USA 99: 123-127. 

Sourjik, V. and Berg, H. C. (2004). Functional in- 
teractions between receptors in bacterial chemo- 
taxis., Nature 428: 437-441. 

Stratonovich, R. (1963). Topics in the Theory of 
Random Noise, Vol. I, Gordon and Breach, New 
York. 

Sykulcv, Y., Joo, M., Vturina, I., Tsomides, T. J. 
and Eisen, H. N. (1996). Evidence that a single 
peptide mhc complex on a target cell can elicit a 
cytolytic t cell response., Immunity 4: 565 - 571. 

Taylor, R. J., Falconnct, D., Niemisto, A., Ram- 
sey, S. A., Prinz, S., Shmulevich, I., Galitski, 
T. and Hansen, C. L. (2009). Dynamic analysis 
of MAPK signaling using a high-throughput mi- 
crofluidic single-cell imaging platform., Proc Natl 
Acad Sci U S A . 

Thattai, M. and van Oudenaarden, A. (2002). At- 
tenuation of noise in ultrasensitive signaling cas- 
cades., Biophys J 82: 2943-2950. 

Tu, Y. and Grinstein, G. (2005). How white noise 
generates power-law switching in bacterial flag- 
ellar motors., Phys Rev Lett 94: 208101. 

Tu, Y., Shimizu, T. S. and Berg, H. C. (2008). Mod- 
eling the chemotactic response of Escherichia coli 
to time- varying stimuli., Proc Natl Acad Sci U S 
A 105: 14855-14860. 



19 



Turner, L., Samuel, A. D., Stern, A. S. and Berg, 
H. C. (1999). Temperature dependence of switch- 
ing of the bacterial flagellar motor by the protein 
ChcY(13DK106YW)., Biophys J 77: 597-603. 

van Albada, S. B., Tanase- Nicola, S. and ten 
Wolde, P. R. (2009). The switching dynamics 
of the bacterial flagellar motor., Mol Syst Biol 
5: 316. 

van Kampen, N. G. (2007). Stochastic processes in 
physics and chemistry., North Holland. 

Wadhams, G. H. and Armitagc, J. P. (2004). Mak- 
ing sense of it all: bacterial chemotaxis., Nat Rev 
Mol Cell Biol 5: 1024 -1037. 

Xing, J., Bai, F., Berry R. and Oster, G. (2006). 
Torque-speed relationship of the bacterial flagel- 
lar motor., Proc Natl Acad Sci U S A 103: 1260- 
1265. 



20 



Supplementary information: 
Noise characteristics of the Escherichia coli rotary motor 

Diana Clausznitzer 1,2 ' 3 and Robert G. Endres 1,2 '* 

1 Division of Molecular Biosciences, Imperial College London, SW7 2AZ London, UK; 

2 Centre for Integrative Systems Biology at Imperial College, Imperial College London, 

SW7 2AZ London, UK; 
3 BioQuant, Universitat Heidelberg, 69120 Heidelberg, Germany 
*e-mail: r.endres@imperial.ac.uk 



Contents 

1 Stochastic differential equations [l| 

2 Model for motor switching 

3 Linearisation of the mo del i 

4 Response functions 

5 Noise spectra Q 

6 Number of high-frequency filters 

7 Alternative Master-equation approach @ 

8 Test of Langevin approximation for motor dynamics @ 

9 Parameters 

1 Stochastic differential equations 

In the main text, we presented a simplified model of the chemotaxis pathway to illustrate signalling and 
noise transmission. Here, we discuss a model for the full signalling pathway and present the response 
functions and noise spectra. Equation (1) in the main text describes the total signalling activity A c of all 
receptor complexes in a cell in response to changes in the methylation level of the complexes and ligand 
concentration is given by 

dA c dA dMj dAdcj 



1 



The dynamics of the methylation level of complex j is described by 



' m ' - 7«(^-^)-7BA^+^ M .(f) (2) 



dt 



lR (N - A,) - ^-A,N 2 Bp + r, Mj (t) (3) 



Note that here we explicitly include the number of CheB-P (B p ) molecules N B in the demethylation 
term, with V cc \\ the cell volume (cf. Eq. (3) in the main text). We denote the activity of complex j by 
Aj. The dynamics of the concentration according to Eq. (2) in the main text is 

dc-j die) , . ... 

In addition, we take into account phosphorylation and dephosphorylation of CheA (A p ), CheY (Y p ) and 
CheB (B p ), which are described by the following equations: 



dN Yp ( k 



dt \ Vcdi 

dN Bp ( k b 



(^-) {N B ,u* - N Bp )N Ap + r) A<P (t) + VA,Y p (t) + VA, Bp (t) (5) 

V ^cell / 

(AV,tot - N Yp )N Ap - k^ y N Yp - r) AtYp (t) + V - Yp (t) (6) 

(A^tot - N Bp )N Ap - k^ b N Bp ~ i lA , Bp (t) + V _ Bp (t) (7) 



dt \ Vcdi 

with Ni the number of molecules of species i in a cell volume. Note, that we neglected the binding 
of CheY-P to its phosphatase CheZ for simplicity, and describe dephosphorylation of CheY-P by the 
effective dephosphorylation rate k- y N Yp . As shown in Sec. [HI this simplification has no qualitative effect 
on the response function. The term r] AtP describes the noise associated with CheA autophosphorylation. 
The terms rj A ^ Bp (t) and T] AtYp {i) represent the noise generated in phosphorylation of CheB and CheY by 
CheA, respectively. The terms r\_ Bp (t) and r]- Yp (t) describe the noise associated with dephosphorylation. 
Note that some noise terms appear in two equations. This is due to the fact that we assign noise terms to 
a specific process, e.g. phosphorylation of CheY by CheA. Hence, the corresponding noise term (t) 
appears in the dynamics of CheA-P and of CheY-P. As a positive fluctuation in the dynamics of CheA-P 
due to phosphorylation of CheY corresponds to a negative fluctuation in the dynamics of CheY-P, these 
noise terms appear with opposite signs in the two equations. The noise intensities and parameter values 
of the model are summarised in Sec. |H1 Finally, the dynamics of the motor is described in terms of the 
probability of the tumbling mode Pew 

- k + (N Yp ) (1 - P CW ) ^ k- (N Yp )P cw + VPcw (t) , (8) 

where we use experimentally derived switching rates k + and fc_ as a function of CheY-P concentration 
(cf. next section and Fig. 2C in the main text). 



2 Model for motor switching 

ITurner et al .] (Il999h presented a model for motor switching to explain the observed motor switching 
rates. The model for motor switching is an MWC model, where 26 subunits of the motor assume one of 
two states corresponding to CW and CCW rotation. While these subunits bind the molecule CheY-P 



2 



independently of each other, the switching of states occurs cooperatively. The authors derive the overall 
rates of switching (averaging over all possible CheY-P occupancy states) as 



k + (Y p ) = fc+(0) • exp i m 



^coop 



■ In 



k-(Y p ) = fc-(O) • exp |m coop 



In 



-f^ccw 
Kqcw 



Y p \ 



Kccw J 
Y 



Kcw 



(9) 
(10) 



where m coop = 26 is the number of motor subunits, — fceTm(/i) is the free-energy difference of switching 
per molecule of CheY-P, and i^ccw and Kcw are the dissociation constants for binding CheY-P in the 
CCW and CW state, respectively. 

Motor switching rates /c+ and fc_ have been derived exper imentally as a funct ion of the concentration 
of a signalling mutant CheY**, which is constitutively active ([Turner et al. ■ Il999h . as shown in Fig. 2C in 
the main text. To obtain the switching rates in terms of CheY-P, rather than signalling mutant CheY**, 
we rescaled the dissociation constants of CheY binding to the mot or such that the sw itching rates are 
equal, i.e. CW bias about 1/2, at CheY-P concentration 3.2 [iM dCluzel et all l200dh. We fitted the 

et all 



above model to the experimental data in Fig. 2C and the CW bias at 33° ([Turner et all 1 19991) and used 
the rate constants k + and fc_ in our full pathway model (Eq. flSJ). 



3 Linearisation of the model 

Similar to the presentation for the simplified model in the main text, we linearise Eq. (HI)-® and insert 
the Fourier transforms of the dynamical variables to obtain the Fourier transformed response functions 
Xr(u>) and noise spectra Sr(uj) for the signalling pathway. The linearised equations read 



d(SA c ) A C1 ^ . dAdlScj) dA 



dt 

d(SN Ap ) 
dt 

d(SN Yp ) 
dt 

d(6Ns p ) 
dt 

djSPcw) 
dt 



dc dt 



dM 



(11) 



X 2 SA C - 


X 3 SN Ap + 


X 4 6N Yp + X 10 SN Bp + r] Ap {t) + T] A<Yp {t) + VA, Bp (t) 


(12) 


X 5 SN Ap 


- A 6 <5iVy p 


- VA,Y p (t) + V-Y„(t) 


(13) 


\uSN Ap 


- XuSNe 


• p - r)A,B p {t) + V-B p (t) 


(14) 


X 7 SN Yp - 


- A 8 <5Pcw 


+ VPcw 0) 


(15) 



with rate constants of the linearised model given in Table [3] 



4 Response functions 

The response functions can be calculated from the linearised Eq. (frT|) - (fT5| without noise after inserting 
the Fourier transforms of the dynamical variables. The Fourier transformed response functions of CheA- 
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P, CheY-P and the motor are 

XN Av (u) : 



dA 

-iLu\ 2 N c — (A 6 - iu) (A12 - iu) 
ac 



(16) 



OA 



dA 

-A4A5 (A12 - iu)] + X 2 X 9 XuN c -^ (A 6 - iu) 



XN Yp (u) 
1p c w M 
AW M 



A, 



A6 — tO,' 

A 7 



As — iu 

An 
A12 - iu 



Xn Av - 



(17) 
(18) 
(19) 
(20) 



From these equations we observe that CheA-P, CheY-P and the motor are in a cascade. At each level of 
the cascade, a new filter proportional to (Aj — icu)^ 1 is introduced which simply multiplies the response 
function of the previous level of the cascade. The characteristic frequencies Aj contain the forward and 
backward rates of the relevant processes. 

5 Noise spectra 

The noise spectra can be calculated from the linearised Eq. f|li p -(|15 |) ■ After inserting the Fourier trans- 
forms of the dynamical variables, calculating the absolute squared value and averaging, we obtain the 
noise spectra for CheA-P, CheY-P and the motor shown in Fig. 3 in the main text 



Sa b (u) = (j(A 3 - iu)(X 6 - iw)(A i2 - iu) - A 4 A 5 (Ai 2 - iu) - Ai Au(A 6 - iu) 



Nr. 



dA 



dc 



u- I — I S c (u) + u z {S a (u) + Q M ) 



+ 



+ 



dA 



^ Nc ~dM ^ 3 ~~ iuj ^^ 6 ~ iujS) ~ A|Ar - 



?-b„ + 



dA 
dM 



Q- 



Q 



dA 

X9XiiN c -^(X e - iu) Qa p + X 4 X 9 XuN c 
dA 

dA 2 \ 

AgAiiA^-— [(A 6 - iu) - A 4 ] Qa,y p I • 



dA \ 

^ 1 'dM t^ 3 - * w )(^6 - iu)(Xi2 - iu) - Ai An(A 6 - iu)- 

dA 

-A 4 A 5 (Ai2 - iu)] + A 2 A 9 AiiA c ^j(A 6 - iu) 



(21) 
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Sn Ap (oj) = ^N c \ 2 2 io 2 [^j \(\ 6 -ilo)(\ 12 -iu J )\ 2 S c (u)+ 

+ |A 2 (A 6 - iw)(Ai2 - iuj)\ 2 N c (u 2 S a (uj) + Q M ) + 

f dA \ 2 
+ ( Ai^: -iwj (A 6 -iw)(Ai2 -iw) Qa v + 

( OA ( c 

+ 



dA ( dA 

(A 6 - iuj) ( -\ 2 \ 9 N C — + Aio ( Ai— - zu; 



Q — B v + 



dA 



+ 



(A 6 - iw)(Ai2 - zw)(Ai— - iw) + (A 6 - iu})(X 2 X 9 N c — + 



Qa,b p + 



dA \ ( dA \ 

dA \ 

Ai - iu j ■ [(A 3 - iuj) (A 6 - iw) (A12 - iuj) - \ w \n (A 6 - iuj) + 

dA 

-A4A5 (A i2 - ioj)} + A 2 A 9 AniV c ^: (A 6 - iu] 



(22) 
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S N (u) 



SpcwM = 



dA 
~~dc 

2 



■ N c X 2 2 u 2 I ^ ) |(A 6 - iw) (A12 - iw)r 



|A 6 — 

+ |A 2 (A 6 - iw)(Ai2 - iw)| 2 iV c (w 2 5 a (w) + Q M ) + 



+ 



+ 



dA 



Xl ~dM ~ iUJ ' ^ e ~ iuJ ^ Xl2 ~ iuJ ^ 



Qa v + 



dA ( dA 

(A 6 - iu) [ -\i>^N c — + Aio ( Ai— - 



A 5 



Ak — iio 



dA 



+ 1 



<2-B P + 



-Y„ + 



+ 



dA dA 
(A 6 - iw)(Ai2 - «w)(Ai— - iw) + (A 6 - zw)(A 2 A 9 iV c — + 



Qa,b p + 



+ 



\6 — ZCJ 



<9A 



- 1 



Q 



A,Y 



dA 



Xl d~M _ / ' ^ Xs ~ ^ ^ 6 ~~ *~ Al2 ~~ ~~ AioAn (A 6 - iu) + 

dA 

-A4A5 (A12 - iu;)] + A 2 A 9 AiiiV c ^: (A 6 - iu) 



X 7^N Yv ( w ) + Qp C v 



I A 



8 - M 



(23) 
(24) 



6 Number of high-frequency filters 

In Fig. 2A in the main text it is apparent that our model does not fully reproduce the high-frequency 
response. The high-frequency response seems to be a third-order filter in the frequency range shown, 
while our model only produces a second-order filter due to CheY-P and motor-switching dynamics. A 
third filter due to the autophosphorylation dynamics, which is included in our model, becomes only 
relevant at higher frequencies. Here, we discuss where an additional filter could originate. 

We explicitly consider the CheY-P/CheZ binding step, and write down the equations for the dynamics 
of the concentration of CheY-P, denoted by y, and of CheY-P/CheZ complex, yz, 



dy 

= 9y (Y tot -y)a- g x (Z tot -yz)y + g 2 yz 



d{yz) 

dt 



= gi {Ztot - yz)y - (g 2 + 33) yz, 



(25) 
(26) 



where a is the concentration of phosphorylated CheA and gi are the rates of phosphorylation of CheY 
(gy), CheY-P/CheZ complex formation (gi), dissociation of CheY-P/CheZ complexes (g 2 ) and CheY-P 
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dephosphorylation (33). Linearising around the steady state (a* ,y* ,yz*) yields 



= g Y (Y tot -y*)Aa~[g Y a*+g 1 {Z tot -yz*)}Ay+{g 1 y*+g 2 )Ayz (27) 
at \ 

Ai A2 A3 

d(Ayz) 



dt 



[giy* + <? 2 + 53] Ayz + gi (Ztot - yz*) Ay. (28) 



A3 +93 A 4 

Hence, we obtain for the Fourier transform of deviations in the CheY-P concentration 

Ag = A|i '-" ^-'l- . . A q. (29) 

(-iw + A 2 )(-iw + A 3 + .93) - A3A4 

To make the analysis easier, we can factorise the polynomial in the denominator, 

a- Xi(-iuj + X 3 + g 3 ) 
(— tu> + ai)(—iw + a 2 ) 



with ai j2 = (A 3 + g 3 + A 2 )/2 ± ^ (A 3 + g 3 + A 2 ) 2 /4 - A 2 (A 3 + g 3 ) + A 3 A 4 - 

We are interested in the behaviour of the frequency-dependent prefactor. Specifically, we ask if by 
considering the CheY-P / CheZ complex formation we obtain an additional high-frequency filter compared 
to Eq. (|18|) . It is obvious from Eq. (|30[) that under most parameter combinations we obtain 1/lu behaviour 
at high frequencies. Hence, no additional filter is introduced. A special case appears for A3 + k3 3> 
01, a 2 . In this case, a l/ui 2 behaviour is observed for medium frequencies max(ai, a 2 ) > w > A3 + g 3 . 
Hence, additional filter appears. At high frequencies to > A3 +.93, the prefactor has 1/cj behaviour. 
However, analysing the expressions for a\ and a 2 reveals that max(ai, a 2 ) is always greater or equal 
to (A3 + g 3 + A 2 )/2. Therefore, this case does not occur for our dynamics. In conclusion, considering 
CheY-P/CheZ complex formation does not introduce an additional high-frequency filter. 

Other p rocesses in the s i gnalling pathway neglected here are the oligomerisation of CheY-P / CheZ 
complexes ( Eisenbach . .2004; Bl at and Eisenbachj . 1996bl lah and a potential slow release of CheY-P from 
the sensory complex as discussed by Blat et al.l ()l998l ) Oligomerisation of CheY-P / CheZ complexes for 



efficient dephosphorylation is similar to CheY-P/CheZ complex formation considered above, and by a 
similar discussion does not introduce an additional high-frequency filter. However, a delayed release of 
CheY-P represents effectively a step between CheY phosphorylation and motor switching in the signalling 
pathway, and hence could introduce a relevant filter if the process is sufficiently slow. Another possibility 
to explain the steep frequency-dependence of the response function is that the duration of experimental 
pulses was long enough to leave a signature. 



7 Alternative Master-equation approach 

Alternative to the Langevin approach, which assumes small fluctuations, we can write down a Master 
equation for the pathway. Here, we focus on ligand and mcthylation dynamics at the receptor cluster. 
Each state of the pathway is described by the variables [Lj, Mj] at each of the receptor complexes, 
where Lj — cjs 3 is the number of molecules in a small volume at the receptor complex and Mj the 
total methylation level of the receptor complex. Assuming for simplicity only one receptor complex, the 
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Master equation for the probability density p is 
dp(L, M, t) 



Of 



k D (s 6 c )p(L-l,M,t) 

+k D {L + l)p{L + l,M,t) 

+ 7R [1 - A(L, M - l)]p(L, M - 1, t) 

+~, B [A(L, M + l)] 3 p(L, M + 1, t) 

-{k D {s 3 c + L) + k D (s 3 co) + 7b[1 - A(L, M)] + j B [A(L, M)] 3 }p(L, M, t) (31) 

For small noise, using van Kampens VL expansion (jvan Kampenl [2007h for the variances of fluctuations 
in the number of ligand molecules and the methylation level at steady state (jAauino et all 120111 ). we 
obtain 

{SL 2 } = c s 3 (32) 
ISM') = 1 | 7fl (3-2A*)(M) 2 c 

x 1 (3-2A*)j8 A* 2 {l-A*)[k D +~f R {3~2A*)(l~A*)P]p' V ' 

with (3 = 1/2 the free-energy difference due to adding one methyl group (in units of ksT). The first 
term is due to fluctuations in the rate of methylation and demethylation and the second term is due 
to transmitted fluctuations in the activity from ligand noise. This corresponds to the results from the 
Langevin approach. Specifically, for ligand fluctuation we obtain the same variance after integration of 
the power spectrum Eq. (17) in the main text. Furthermore, the power spectrum of the receptor complex 
methylation level using the simplified model from the main text is 

sm{uj) = H^? M • (34) 

The variance of the methylation level, obtained by integration of the power spectrum, corresponds to 
the above result. 



8 Test of Langevin approximation for motor dynamics 

We chose to describe the dynamics of the motor using the Langevin Eq. (4) in the main text 

= fc+(l-P cw )-A : _Pcw + ?7Pcw(i) (35) 

at 

with switching rates from CCW to CW (first term) and from CW to CCW (second term), as well as an 
additive Gaussian white noise term (last term) with zero mean and autocorrelation (?7p cw W^Pcw^')) = 
Qp CVJ 8{t — t') with Qp cw = 2fc + (l — P t *) = 2fc + fc_/(fc + + fc_). For constant switching rate constants k + 
and fc_ , the power spectrum Pew is (cf. Eq. (9) in the main text) 

"""■M- sJCl*-? - (36) 

To see that this is a valid description o f the binary motor-s witching process, we calculate the spectrum 



exactly according to the derivation bv IStratonovichl ( 19631 ). For a stochastic two-state process, whose 



time interval lengths in each of the two states t± and t%, respectively, are independent and identically 
distributed random variables, the power spectrum is given in terms of the Fourier transforms of the 
waiting time distributions 0i(cj) and 02(w) for each of the states, 

q(t , 2 [i-exHHi-^H] 

5M -^(( Tl) + ( T2 )) Jt i- 01 He 2 H ' ( ' 7) 



s 



0.5 



• two state 
x Langevin 




2 4 6 8 10 

Fold Q Y 

Figure 1: Variance of the motor bias as a function of CheY-P noise intensity aQ Y — 2ak Y for Langevin 
and two-state dynamics. 



where 5R indicates the real part. Assuming for the motor that switching between the sta tes CW and CCW , 
respectively, foll ows exponential interval distributions determined by rates fc + and ( Blo ck et al. . 1983t 



Bai et al.L l201oh . the Fourier transforms of the waiting time distributions are given by, 



e C wM = r— j (38) 

fc + — IU1 

&ccw{u) = fc ~ . , (39) 

and the power spectrum is 

S2{uj) = (fc + + M- 2 + (fc + + fc-) 2 ' m 

This result is equivalent to the spectrum obtained from the Langevin equation. Furthermore, we tested 
numerically that the statistics of the Langevin equation and binary process are equivalent for fluctuating 
rate constants k + and fc_ due to the CheY-P dynamics. We simulated time courses of CheY-P according 
to the simplified equation 

dN Y 

-^ = k Y -k_ Y N Yp + VYp (t) (41) 

with rates k Y — 5/s and k Y such that (Ny)* /V cc \\ = (Y p )* = 3.2/iM. The noise term r] Yp (t) is Gaus- 
sian and white with zero mean and autocorrelation (r] Y (t)rj Y (t')) = 2k Y aS(t — if) = Q Y aS(t — <'), 
where we varied a. Fluctuating CheY-P was translated into the rates k + (Y p ) and fc_(Yp) according 
to Fig. 2C in the main tex t. The Langevin equation was simulated using a Euler-Maruyama algo- 
rithm ( Kloeden and Platenl . 1992 ) and the binary process using a Gillespie algorithm. Figure [T] shows 



the variances of both processes as obtained from 10 2 runs for each value of a. As can be seen from the 
figure, the Langevin equation is a good description for the binary process of motor switching. 



9 Parameters 

Rate constants and total cell concentrations of proteins for the full pathway model are given in Table [T] 
The noise terms f]Aj, Tj C s, T)Mj an d ??p cw are the same as in Eq. (15) and (18)-(20) in the main text and 
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their power spectra are given there. The noise associated with phosphorylation and dephosphorylation 
VAp, i]a,b p ,Va,y p , V-b p and ?7-y p are assumed to be Gaussian white noise terms with zero mean and 
autocorrelations (r]i(i)r]i(t')) — QiS(t — t') with noise intensities Qi given in Tabled The linearised rate 
constants for the full pathway model are given in Tabled Fitting parameters of the Fourier transformed 
response function Fig. 2 are listed in Table |4] Parameters for Fig. 4 are listed in Tables [5] and [6l and 
those for Fig. 7 are listed in Table [JJ 



Table 1: Parameters of the full pathway model, including references to literature. The literature values 
are given in parentheses where different from our parameter values. fc_ y was determined by the condition 
that at steady-state with A^=l/3, the concentration [Y p ]* = [V] to t/3 (pouriik and Beret 120021 ). 



Parameter 



[A] iot 

[B] tot 
[Y] toi 

V ce ll 

Na.ua 

N B ,UA 

N Y ,tot 

JVtot = NN C 

k 2 

k,\ 

ky 

kts 

k-Y 

k-B 

1R 
IB 



Value 



5 (M 
0.28 
9.7 fxM 
1.4 fl 
4215 
236 
8177 
7027 
10 3 s" 1 
10 s- 1 

100 fiM' 1 a- 1 
15 /dVI" 1 s- 1 
5 s- 1 
1.35 s^ 1 
0.006 s- 1 
3.14 ^M- 2 s- 1 



Reference 



Souriik and Berg (2002) 
Li and Hazelbaucr (2004 
Li and Hazelbaucr (2004 
Souriik and Be rg (2002) 
calculated from above 
calculated from above 
calculated from above 
Souriik and Berg (2002) : 
Shapovalov and Lester 
Wolanin et al. (2006) 
Stewart et al. (2000) 
Stewart et al. (2000) 
adjusted to yield steady-state value 
(0.35 s" 1 ) Brav and Bourret (1995) : Stewart 
Clauszni tzer et al. (2010) 
Clausznit zer et al. (2010) 



; Li and Hazelbaucr (2004) 
2004) 



Table 2: Intensities of Gaussian white noise terms in the full pathway model. Index i represents noise 
term rji. 



process 


index i 


noise intensity Qi 


receptor switching 


a 


2k 2 A* 


ligand diffusion 


L 


2Dsc 


receptor de/methylation 


M 


2 1R {N-A*) 


CheA autophosphorylation 


Ap 


K (itv) (N A , tot - N Ap ) 


CheY phosphorylation 


A,Y P 


Gfe) (Ny,tot-N Yp )NX p 


CheB phosphorylation 


A, B p 


fe) (NB,tot-N* Bp )NX p 


CheY dephosphorylation 


-Y P 


k-yN Yp 


CheB dephosphorylation 


—Bp 


k- b N* Bp 


motor switching 


Pew 


2k* + k*_ 
k* + +k'_ 
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Table 3: Parameters of the linearised equations for the full pathway. 



A, 



expression 



Ai 
A 2 

A 3 

A 4 

A 5 

A 6 

A 7 

A 8 

A 9 

Aio 

An 

Al2 



IR 



N1J 



a: 



(l&r) + (fe) (^.tot ~ Np f ) + (^) (N B , tot - N* Bp ) 



k Y N 

Vcoll , 


NX 


k Y ^ 

Koii / 


(JVy.tat - 


ky ^ 
Koll , 





i / /-i _ p*\ dk + _ p* dk- 
V co „ r ' dY v r dY v 







VoeU 






| fee 

V Vcell / 


AT* 

yip 




( k B 
V Veil , 


(iVs.tot - 


N* Bp ) 


| fee 


yip 


B ■ 



Table 4: Fitting parameters for response function of the full pathway m odel for Fig. 2 in th e main 
text. Motor switching rates where not adjusted when fitting to the data bv lShimizu et al. ( 2010h as the 
high-frequency response was not measured in these experiments. 



Parameter 


Block ct al., 1982; Scgall ct al., 1986 


Shimizu et al., 2010 






32° C [s" 1 ] 


22° C [s- 1 ] 


adaptation: 








Xi(dA/dM) 


0.178 


0.018 


0.0039 


A 9 


0.0263 


0.0027 


5.6 10- 4 


motor switching: 








A 7 


4.4 10" 4 






As 


2.111 







Table 5: Parameters for cell-to-cell variation in Fig. 5 A in the main text. For additional parameters sec 
Tabled! 



Parameter 


WT1 

(black line) 


red line 


green line 


blue line 




1.05 


52.4 


1.05 


1.05 


k*_ [s- 1 ] 


1.06 


53.0 


1.06 


1.06 


IR [S X ] 


0.0069 


0.0069 


6.9 10" 5 


0.0069 


7s [MM" 2 s" 1 ] 


3.14 


3.14 


3.14 10^2 


3.14 


N to t 


7000 


7000 


7000 


70 
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Table 6: Parameters for cells with low motor bias in Fig. 5B in the main text. For additional parameters 
see Table [T] S imilar spectra wit h increased low- frequency component have been found for cells with other 
motor biases ( Park et al. . 2010h . 



Parameter 


value 


N tot 

Y tot \jjM] 
k* + [s _1 ] 

dk + /dY p [s- 1 ^M- 1 ] 


700 
2 

0.015 
4.75 



Table 7: Parameters for the fluctuation-response theorem in Fig. 7 in the main text. For additional 
parameters see Table [TJ 



Parameter 


varying j R (solid) 


varying -jb (dashed) 


varying [Y) tot (dotted) 


iV to t 


700 


700 


700 


1R IS" 1 ] 


varied 


0.0069 


0.0069 


IB [S- 1 ] 


3.14 


varied 


3.14 


[*1tot [mm] 


9.7 


9.7 


varied 
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